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INTRODUCTION 


This subject program was instituted in FY 1965. It is a cooperative endeavor of the 
American Society for Engineering Education, the University of Maryland, the Catholic 
University of America, and Goddard. This program is intended to operate on a two-year 
cycle. Professors agree to participate during two consecutive summers. 

The objectives of this program are many; some of them are listed below: 

1. Stimulation of schools to become interested in the research problems now con- 
fronting Goddard. 

2. Creation of interest on the part of professors to continue their research after 
completing the formal program. 

3. Stimulation of our people professionally through associations with the professors 
and through their participation in the program's seminars. 

4. Establishment of closer ties wittt the Universities. We have been able, through 
this program, to develop working relationships with the University of Maryland 
and Catholic University which are outstanding. This is important for we have 
more than 150 employees attending these schools. 

5. An advertising tool for Goddard which is most acceptable professionally. Every 
major educational organization has been exposed to the fact of our participation. 
ASEE is an influential professional organization. 

Now that we have completed the first summer of this program we are happy to re- 
port that the progress made on the above objectives was beyond our expectations con- 
sidering that this was the first summer. 

One of the goals that has been implied is the solution or approach to the solution of 
Goddard problems that might be of interest to the Summer Fellowship participants: The 
problems that they selected are listed in the table of contents and the results are com- 
piled in this report. It should be borne in mind that the time allotted for the solution of 
the problem was very short (an equivalent of seven weeks) and consequently some of the 
problems need more time for complete solution. Most of the 1965 Summer Fellows are 
expected to return for the summer of 1966 in addition to the 15 new Fellows that we hope 
to have this summer. 



John T. Mengel 
Assistant Director, Tracking & Data Systems 
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v R DESIGN 
William D. Stanley 

Assistant Professor, Clemson University 
NASA-ASEE Fellow 


I. INTRODUCTION 

Large-scale processing of telemetry data frequently requires filtering processes to 
remove undesirable signals and noise that perturb the signal. In the case of analog sig- 
nals, this process is normally accomplished by the use of RLC passive or RC active fil- 
ters whose effect on the signal is usually described by the Fourier transform of the im- 
pulse response. This characteristic is called the frequency response and consists of an 
amplitude response and a phase response. 

In processing sampled or digital data, it is often inconvenient or undesirable to em- 
ploy analog filters for such filtering operations. Rather, it seems more feasible to pro- 
gram a numerical algorithm on a digital computer which accomplishes the same task on 
the digital data. The computer then becomes an integral part of the overall system, and 
furthermore, it has the extra advantages of flexibility and ease of adjustment and change 
that is usually not present in an analog filter. 

The purpose of this paper is to present a study of the design of digital filters. Both 
the approximation and the realization problems will be considered, and techniques will 
be presented for designing both low-pass and bandpass filters. A survey of the literature 
dealing with the fundamental theory will be made, and this will be followed by step-by- 
step design procedures. The resulting "design" in each case .will be a numerical algorithm 
whose effect on the data is the desired filtering operation. 


H. BASIC THEORY 


For the purposes of this paper, the term "digital filter" will be used to designate any 
numerical operation that can be performed on a set of data to alter the spectral repre- 
sentation of the data in some desirable manner. Theoretically, such a filtering operation 
could be done by pencil and paper or by a desk calculator if the proper sequence of oper- 
ations is known. A more practical approach is to program the filter algorithm on a digital 
computer and allow the computer to perform the filtering operation. # The filter operation 
can be performed on-line in real time or the data may be stored on tape and processed at 
a lower rate. ; 

i 

The underlying approach to any digital filter is the theory of sampled-data systems. 

There is an extensive amount of literature pertaining to the subject, and no exhaustive ; 

survey of this literature will be attempted. Among the texts pertaining to sampled-data > 

systems are Ragazzini and Franklin (1), Jury (2), Tou (3), Monroe (4), and Kuo (5). All j 

of these books contain extensive bibliographies of the research literature pertaining to ■ 

sampled-data systems. 1 

This paper will draw basic theorems and results freely from the above-named refer- * 

ences. The proofs of these results will be repeated only when necessary for clarity or ' 

further development. ! 
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A sampled signal may be represented*!* the product of the continuous signal and a 
periodic pulse train. As far as a numerical process is concerned, the pulse train may 
be considered as an impulse train. Let x(t> represent a continuous signal and x*(t) 
represent the sampled signal. The sampled signal may be represented as 


x*(t) - x(t)i T (t) 

where t is the time between samples and 


(1) 


co 

22 & ( t - nT ) 


( 2 ) 


\ 


as illustrated in Figure 1. 


x(t) 
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x*(t) 


ill 




Figure 1 


Equations (1) and (2) are equivalent to 

m 

x » (t ) = 22 x < nT ) h - nT ) 


(3) 


Let f represent the process of Laplace transformation. The following definitions 
are made: 


C fx(t>] = X(s) 

(4) 

C[x*(t)] = X*(s) 

(5) 


2 


l 


Transformation of Eq. (3) yields 


X*(s) -- 


00 

V I 

/ *("T) 


It can also be shown that (6) 4 s equivalent to 


X*(s) -- - 
T 


L 



( 6 ) 


(7) 


According to Eq. (7) the spectral representation of X*(s) is periodic. This nsa is 
illustrated in Figure 2 for s = jo. (Fourier transform). As long as the minimum Nyquist 
sampling rate is maintained, the original spectrum of x (t) is preserved, although modi- 
fied by a factor, and the same spectrum is repeated on both sides at multiples of l/T. 


ZVA 


X(jc 


UJ 



x*(jw) 




Figure 2 
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The z-transform is ordinarily introduced to simplify the operation of a numerical 
process. Let 




Z r t T * 

(8) 

and 

X(z) i 

(X*(t)! : X(S) ] T# 

(9) 

thus 

X(z) -• 

Y n x(nT) z* n 

(10) 


o 
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Now assume that x*(t> (or equivalently X(z) ) is applied to the input of a linear 
numerical processor as shown in Figure 3. Let y*(t) represent the output and iet 

Y(z) :• -i [y*( t) ] (11) 

It can be shown that 

Y(z) = G(z) X(z) (12) 

Where o<z) represents a type of transfer function. In the case where the processor is a 
hardware item such as an electrical network, G(z) is called the pulse-transfer function. 

In the case where the processor is a software item such as a computer program (the case 
of primary interest in this paper), G(z) is simply an algorithm whose form will be dis- 
cussed in a later section. 


x(i) 


X 'x*(tT 


G(z) 




y*(t) 


Figure 3 


The transformation of Eq. (8) maps the left-hand half of the s-plane to t»n inside of 
the unit circle in the z-plane and the right-hand half of the s-plane to the otunide of the 
unit circle in the z-plane. In the sam 'nner as the spectrum of x*(s) as predicted by 
Eq. (7) is periodic, so is the spectrum of y *(s) periodic. Thus 

v<» -J £ v (• . i ) (13) 

This periodic property is an inherent characteristic of the sampling iro-ess and must be 
clearly interpreted in any numerical process. Basically it is relincu o the familiar un- 
certainty principle. Since there is uncertainty in the time function between samples, 
there is a corresponding ambiguity in the frequency domain representation. 

All this means that it is impossible in a sense to create any digital filter with an 
arbitrary form of characteristic at all frequencies. Instead, the approximation must be 
done in a "local" sense with the expectation that the given frequency function will repeat 
itself at other frequencies in its present sense. This does not mean that the digital filter 
is more limited in its utility than the analog filter. On the contrary, a sufficiently high 
sampling rate may be chosen along with accurate digital to analog conversion in such a 
manner as to essentially eliminate the unwanted portions of the spectrum. Furthermore, 
complex filters may be synthesized quite readily on a digital computer m cases where the 
corresponding analog filters would be very difficult to build. 


III. FILTER REALIZATION FROM TRANSFER FUNCTIONS 

It is normally advantageous to study the realization problem before considering the 
approximation problem, although an actual design usually proceeds in the referse order. 
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It cat. be shown that any physically realizable linear numerical transfer function can 
be written as the ratio of two polynomials in z~ k of the form 


G(z) ^ 


Y(z> 
X(z) ' 


N 


k*0 


«k 2 


-k 



(14) 


For a stable system, the poles of G(z) must lie within the unit circle in the z-plane. 
A more intuitive form of Eq. (14) can be obtained by writing it in the form: 


II 


“ N 

1 * 7~! ^ z ~ k 

•= X(z) 

L 

V-Y 


k 0 


(15) 


The quantity z* k operating on any F(z) corresponds : d delaying f(t) by kTsec. More 
generally 


iz* k F(z)l f(nT-kT) 

Performing the inverse z-transformation of Eq. (15) yields 



M 

y<nT> * L^ y(nT - kT) 

N 

T \ x(n T - kT) 


k-1 

k,0 

or 




N 

>1 


y(nT) 'y"' x(nT - kT) - 

k 0 

J2 l\y(nT-kT) 


(16) 


(17) 


(18) 


Thus, the transfer function of Eq. (14) is equivalent to the numerical operation of Eq. (18). 
The reason for the unity term in Eq. (14) is now clear. A block diagram representing the 
sequence of operations expressed by Eq. (18) is shown in Figure 4. 



Figvf# 4 
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The form of synthesis (or programming) expressed by Eq. (18) is called a "direct 
progr amming " scheme. There are a number of other techniques which are discussed in 
the texts referenced earlier, along with the direct method, but they will not be discussed 
here, ft is felt that the direct method is the easiest technique to implement directly on 
the computer, although some of the other techniques may yield fewer computer opera- 
tions in some instances. 

As an illustration of this method consider the transfer function given by 


G(z) = 22 * 5 

8z 2 + 6z + 1 


(19) 


The first step in realizing G(z) is to put it in the form of Eq. (14). Dividing numerator 
and denominator by the value 8 and multiplying numerator and denominator by z' 2 re- 
sults in 


G(z) - 



( 20 ) 


Thus from Eq. (18), 


y(nT) = 2 x(nT) + - x(nT - T) + 1 (nT - 2T) 

8 ^ 8 

- 2 y(nT - T) - I y(nT - 2T) (21) 

4 o 

ft desired, the notation may be modified to read 

y(n) = -2 x(n) tix(n-l) +2 x(n - 2) 


--f- y ( n “ “4~ y ( n “ (22) 

4 8 

In subsequent work, the notation of Eq. (22) will be employed with the sampling time T 
implied. A flow diagram for Eq. (22) is shown in Figure 5. 


IV. ANALOG TO DIGITAL FILTER TRANSFORMATION 

Now that the basic form of the numerical transfer function has been stated, the next 
question is: How does one relate the basic specifications of a desired filtering operation 
to the final transfer function to be realized? This problem, of course, is the familiar 
approximation problem, and the solution is evidently not unique. One possible approach 
would be to attempt to approximate the transfer function directly in terms of a ratio of 
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Figure 5 


polynomials in z _1 (or equivalent e- sT ), However, this approach appears to have re- 
ceived very little attention due to the difficulty of working with exponential functions. 

On the other hand, there is a wealth of literature available on the approximation 
problem of analog filters. K one can obtain some direct relationship between the analog 
domain and the digital domain, the approximation problem can be carried out in the ana- 
log domain, and then a suitable transformation may map the function into the digital 
domain, thus achieving the desired result. 

The most intuitive approach to this problem is to find an analog transfer function 
representing the desired response and then compute the z-transform corresponding to 
this function. However, as pointed out by Steiglitz (6), (7) and Hauptschein (8), this 
method yields rather poor results due to poor correspondence between a sampled- input 
and a continuous input applied to the same structure. Consequently, this method will not 
be pursued any further. 

It appears that the use of the bilinear transformation is the most fruitful approach to 
correlating the analog domain to the digital domain. The use of this transformation in 
studying digital filters and control systems has been investigated by Johnson, Lindorff, 
and Nording (9), Lewis (10), Carney (11), and all the references previously mentioned in 
this paper (1-8). 

The basic form of the bilinear transformation is given by 


or the inverse transformation 


p = K 


(« - 1) 
z + 1 


z=LL£ 

K -p 


where K is a normalizing constant. 


(23) 

(24) 


f 
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It can be shown that the left-hand half of the p-plane is mapped into the interior of the 
unit circle in the z-plane; whereas, the right-hand half of the p-plane is mapped to the 
outside of the unit circle in the z-plane. In other words, there is a similar correspond- 
ence between the p-plane and z-plane as is true with the s-plane and z-plane. On the 
other hand the relationship between the p plane and z-plane is a simple ratio of first- 
degree polynomials and is not complicated by exponentials. 

The basic procedure for using the p-plane approach is to carry out the basic approx- 
imation in the p-plane using the conventional analog approximation technique. Once the 
desired function of p has been obtained, one merely makes the substitution of Eq. (23) and 
manipulates the expression into the final desired form of a numerical transfer function. 

Normally the approximation process is performed along the imaginary axis in the 
s-plane. It can be shown that the imaginary axis of the artificial p-plane maps to the 
imaginary axis of the s-plane, and thus the normal approximation will usually be per- 
formed along the imaginary axis of the p-plane. An interesting point of consideration is 
the relationship between the imaginary axes of the two planes. In this context, let 


s = a + jw = transform variable of final interest (25) 

p = a + j \ = transform variable in artificial p-plane (26) 


since z = e sT , from Eq. (23) we have 


P - K 


e sT - 1 
e sT + 1 


= K tanh 


sT 

2 


and 


s tanh'" 1 — 
T K 


(27) 


(28) 


For s = j w, p is imaginary and thus 


X - K tan ~ (29) 

and 

w = il. tan* 1 Ji (30) 

T K 

The transformation between x and w of Equations (29) and (30) is illustrated in Figure 6. 
The periodic nature of the resulting transfer function can be deduced from the multivalued 
nature of Eq. (30). Any form of response in x will be repeated at certain multiples of the 
sampling rate along the axes. 

In any particular approximation problem the constant K is determined from either 
Eq. (29) or (30) such as to map any particular x to correspond to any particular w. For 
example, the given analog function of p might be normalized for a cutoff frequency of 
1 rad/sec. Then for a desired cutoff in w, Eq. (29) is solved for K. Thus, if Xr is to 
correspond to wr ; K is determined from 
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Figure 6 


K = 


\r 


tan 


w T 


(31) 


The bilinear transformation has been shown to possess another interesting feature 
that will be considered in the next section. ft will be shown that if the trapezoidal inte- 
gration rule is used to approximate an analog system, the resulting numerical operation 
will be identical with the bilinear transformation applied to the analog function. The 
reader interested only in practical filter design may omit the next section without loss 
of continuity. 


V. RELATIONSHIP OF BILINEAR TRANSFORMATION TO 
NUMERICAL INTEGRATION 

In this section, the following theorem will be developed and discussed: 
Theorem 

The bilinear transformation from the p-domain to the z-domain of the form 


p = 1 L=A (32) 

T z + 1 ! 

I 

is equivalent to applying the trapezoidal integration rule to the set of n-simultaneous dif- 
ferential equations that describe the corresponding analog system. 

In order to prove this theorem a few techniques will first be discussed. The first is 
the well-known trapezoidal rule of integration. Consider a single first-order equation of 
the following form: 


I 


9 



f 


dy « dx 


(33) 


The best approximation in the first-order sense to Eq. (33) is achieved from the equation 


y n -y a -l b (yn + y n .l> ( X n" X n -l> ( X n+ X n- 1 > 


+ a. 


(34) 


where T is the sampling period. Eq. (34) is a general form of the trapezoidal rule for 
the first-order linear differential equation. The solution for y proceed 


or 


now let 




a oM / a 0 T 


y n 


H) s Ml. 

HP K)‘ K) 



n - 1 


(35) 


(36) 


Y(p) = C [y(t>] 

(37) 

X(p) =£ [x(t>] 

(38) 

Y(z) = 9 . [y(n)] = ? [y n ] 

(39) 

X(z) = -i [x(n)] = > [xj 

(40) 


Assuming that the system is initially relaxed, application of the Laplace transforma- 
tion to Eq. (33) yields 

PY(P) + bY(p) = p X(p) + a 0 X(p) 

The transfer function g(p) is 

Y(p) _ a iP + a o 


G(p) = - 

X(p) p + b 

With the same initial assumption, application of the z-transformation to (35) yields 

w [(* * ")- (‘ - ?)*”] = «*{(”' * T-)' (“■ ' ~ 2 -} [ 


(41) 


(42) 
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The transfer function H(z) is 


H(z) = 



(43) 


The transformation between p and z can be determined by equating the righthand sides 
of Equations (41) and (43). Thus 


a iP + a Q 

p + b 



] 

(44) ! 


After some manipulation, the following result is obtained: 


P=i-fLli 
T 2 + 1 


(45) 


which is clearly the bilinear transformation with K = 2 /T# With this choice of K , the 
relationship between x and w is 


A = 



wT 

2 


(46) 


| 

i 


t 

I 

j 


For very low frequencies, 


u 

X £ w for w < < — 
T 


i 

(47) | 


Thus, the trapezoidal approximation is seen to yield an approximate one-to-one cor- 
respondence between the analog domain and the digital domain for very low frequencies. 

The quantities H(z) and G(p) are related by 


h<i,=g (4Ht) 

as proposed. 


(48) j 
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The theorem originally stated has thus been proved for the first order case. The 
extension to the more general order can best be achieved by means of the state-variable 
approach. The state variable formulation represents a system of order n by n first 
order equations of the given form. It can be shown that one possible form for expressing 
such a system can be achieved by the techniques to follow. 


It can be shown that a lumped, linear, time-invariant system can be described by a » 

differential equation of the form. 


d"y cF'V 

+t> n-l 7 + 

dt" dt "' 1 


• V = 


d°* 

dt" 


+ a. 


n- 1 


d"-*x 

dt "* 1 


+ 


+ a o x 


( 49 ) 


which is equivalent to the transfer function 


a n p" + a p"' 1 + 

G(p) = -2 1 i 

P* +b n-i P"' 1 + ' 



( 50 ) 


Let 


d"~*y 

dt "* 1 


d"' 2 y 

dt ”' 2 


( 51 ) 


y„ = y 


The following constants are defined: 


A n = a n 


\-i " B n-1 “ 

K-2 = a n-2 “ ^n-1 

Aj = a l " b n .I Aj “ b n-2 A 3 K 


( 52 ) 


r 


* 
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The system can then be described by the equations 


^1 

dt 


+ D 


n~l 


y t 


* b n-9 ^ 


2 J 2 



+ 


X 



dx 

dt" 


( 53 ) 


dt 


y„-j = ^ 


dx 

IFF 


Let 




y = 



(54) 



n-1 


0 0 


[B] = 


0 



(55) 


0 


0-1 0 



Then, the set of equations of (53) are equivalent to 


dy 

dt 


4 


[B] y = A^ + A 0 x 


(57) 


which places the entire system in matrix form in a manner analogous to that of the simple 
first* order equation. The extension of the theorem can be done on the individual equations 
as in (53) or on the matrix form in (57). 
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VI. PRACTICAL "LOW-PASS" DIGITAL FILTER DESIGN 

The problem of the practical design of low-pass digital filters will be considered in 
this section. Strictly speaking, such a filter should actually be termed a pseudo- low-pass 
filter due to its periodic nature as previously discussed. The form of a typical charac- 
teristic is shown in Figure 7. Thus, if a true low-pass action is desired, the sampling 
rate must be considerably higher than the "cutoff" frequency, and accurate digital-to- 
analog conversion must be employed. 



Using the bilinear transformation, the technique for synthesizing a filter from the 
transfer function can be summarized as follows: 

1. Using the theory of the analog approximation problem, obtain a transfer function 
for the low-pass filter of the desired form in the analog domain. In this context, 
the ITT Handbook (12) contains response curves for the Butterworth, Chebyshev, 
and maximally flat-time delay filters, and the text by Weinberg (13) contains 
normalized element values. Any of the other texts on network synthesis can be 
employed in this task. 

Thus, by some standard approximation method, a function of p , say G(p) , is 
obtained as a ratio of polynomial representing the desired response in the analog 
domain. It is not necessary (or desirable) to denormalize G(p) at this point such 
that the cutoff frequency in the analog domain has some unique correspondence 
to the cutoff frequency in the digital domain. Rather it is simpler to develop a 
normalized G(p) such that p = j a. = j l is the reference frequency in the analog 
domain. The total process of normalization results from step 2 to follow. 

2. The second step is the determination of the normalization constant k. Assume 
that some particular reference analog frequency, \ r , is to correspond to some 
final real reference frequency wr. The quantity K is determined from the 
relationship 

K = Ar cot (58) 

3. With k known, the final transfer function h(z) is determined by substituting 


p = K 


z - 1 

Z + 1 


(59) 
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into G(p) . In other words 


H(z) = G(p) (60) 

(* + l) 

4. The quantity H(z) is arranged in descending powers of z and synthesized accord- 
ing to the technique of Sec. m. 

5. If a plot of the frequency response of the digital filter is desired, the transforma- 
tion curve of Figure 8 can be used as will be discussed in the next section. Notice 
that the abscissa of the curve is the normalized quantity f/f T , where 


2T 


= haif of sampling rate 


(61) 



Flgur* 8 
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vn. LOW PASS DESIGN EXAMPLES 


Example 1. A second-order filter of the Butterworth or maximally-flat class will 
be synthesized with the following specifications: 

Sampling Frequency = f s * 1 'T = 2 kc 

Cutoff Frequency - f c = 500 cps = 0.5 kc 

The analog response for such a filter with kc- 1 rad/sec is shown by Van Valkenburg 
(14) to be 


G(p) =— J 

P 2 + v' 2 p + 1 


(62) 


The constant K must now be determined such that \ = 1 maps into »«2»x 500 rad/sec. 
Thus, from Eq. (58) 


K cot ~ S £2ZL = j (63) 

2000 

The desired transformation is 

P = ^ (64) 

2 + 1 


which substituted into Eq. (62) yields 


H(z) 


Z 2 + 2 Z 4 1 

(2 t v 2 ) z 2 + (2 - /2) 


(65) 


This can be puMn the proper form for programming by dividing numerator and denomi- 
nator by (2+ v'2)z 2 = 3.4142. Thus, 


H(z) - 0-29289 (1 + 2*’ 1 ± «' 2 ) 
1 + 0.82841 z* 2 


( 66 ) 


The ajqiropriate algorithm is 


y„ = 0.29289 (x n + 2x n ., + x 0 . 2 ) - 0.82841 y n . a (67) 

A curve of the frequency response may be plotted by using the ITT Handbook curves and 
the transformation curve of Figure 8. A number of points are tabulated in the following 
table. The quantity o<db) represents the input-output attenuation in decibels obtained 
from the curves, and the quantity G represents the magnitude of output to input. A plot 
of several "cycles" of the response is shown in Figure 9. In this case the periodic 
nature of the response is quite apparent due to the low-ratio of sampling frequency to 
cutoff frequency. 
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A *3 

.4 

.5 

.6 

.7 

.8 

.9 

1.0 

1.5 

2 

2.5 

ii 

G(db) .05 

.11 

.28 

.54 

.94 

1.48 

2.15 

3 

8 

12 

16 

i 

G 


.968 

.94 

.897 

.84 

.78 

.707 

.395 

.25 

.158 

0 

f/lkc 


.3 

.345 

.39 

.43 

.47 

,5 

.623 

.705 

.76 

1.0 


Note: Table is accurate only to slide rule and curve accuracy. 



Figure 9 


Example 2. As a more complex example, a fourth-order Chebyshev filter will be 
synthesized to meet the following specifications: 

Pass-band Ripple = « = 1 db 

f, = 2 kc 

f e = 200 cps = 0.2 kc 

The "cutoff’ frequency in this case is considered to be the last 1 db-down frequency. 

For w c = 1, the analog transfer function is 

G(p) = i <68) 

1 + 2.6943 p + 5.2748 p a + 3.4568 p 3 + 3.6280 p 4 

The constant k is determined from Eq. 58 with and *r = 2n x 200 

K - C ot = cot lo° * 3.0777 (69} 
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The transformation is 


p = 3.0777 

z + 1 


(70) 


Substitution of Eq. (70) into Eq. (68) yields after a lengthy but straight-forward 
computation 


H(z> = _ 
1 


0.024717 (1 + 4z -1 + 6 z~ 2 + 4 z~ 3 + z~ 4 ) 
3.0S43z' 1 + 3.8290z" 2 - 2.2925z -3 + 0.550742" 4 


The associated algorithm is 


(71) 


y n = 0.024717 


[x n + 4x n-l + 6x „-2 + 4x „-3 + X „-4 ] 


+ 3.0543 y n _j 


- 3.8290 y n . 2 + 2.2925 y n _ 3 - 0.55v74 y n _ 4 (72) 

A table consisting of a few points of the frequency response is shown below, and a plot of 
this response is shown in Figure 10. 


K 

K 

0 

.383 

.707 

.924 

1 

1.05 

1.58 

2.1 

2.63 

a> 

*3.0777 

0 

.124 

.23 

.3 

.325 

.341 

.513 

.683 

.855 

00 

f/lkc 

0 

.08 

.14 

.18 

.2 

.21 

.305 

.38 

.455 

.5 

G(db) 

1 

0 

1 

0 

1 

3 

24 

36 

45 

00 

G 

.891 

1 

.891 

1 

.891 

.707 

.063 

.016 

.006 

0 


Note: Table is accurate only to slide rule and curve accuracy. 


*$ince p = 3.0777 


z- 1 
z + 1 


= 3.0777 torch y, 


A. 


= 3.0777 tan 


wT 

2 


or 


3.0777 


= ton 


wT 

2 


VIE. BAND-PASS FILTER TRANSFORMATION 

The techniques under consideration may be used to design a pseudo-band-pass filter 
of which a typical frequency response curve is sketched in Figure 11. As "in the low-pass 
case, the response is periodic along the w-axis. The point that characterizes the band- 
pass case is the fact that in the first "cycle," the pass-band occurs somewhere removed 
from the lowest frequencies. 

One way to approach the design of a band-pass filter is to first obtain the transfer 
function of an analog band-pass filter and employ the method of the preceding two sec- 
tions. However, this method may require a lengthy set of computations in the analog 
domain, and the transformation in maintaining correct bandwidth is more difficult to 
interpret. Rather a method will be developed that will yield directly a band-pass digital 
filter from a low-pass analog prototype by means of an extension of the basic bilinear 
transformation. 




f/lkc 
Figure 10 


" AMPLITUDE RESPONSE 

a. uximuL . . 

1/2T 1/T 3/2T 2/T 

Figure 11 
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Assume then that a given low-pass prototype is known in the analog domain with 
transfer function G ( p ). It will be shown that the following transformation will yield a 
pseudo-band-pass digital filter. 


P = A 




z - 1 
z + 1 


+ 


a 2 (z + 1) 
(z - 1) 


(73) 


= A 


. . ST 2 . ST 

tanh — + a* coth 

2 2 


(74) 


The relationship between the w-axis and the A. -axis is obtained by letting s = j w and p = j k. 
The result is 


K = A 




(75) 


W 


Let (76) 

a 

Then Eq. 75 reduces to \ = aA j^x j (77) 

which is equivalent to the normalized low-pass to band-pass transformation of analog 
filter theory. Sketches of Equations (75) and (76) are shown in Figure 12. 




Figure 12 


From Figure 12b, it can be observed that for a low-pass reference frequency \r and 
the corresponding negative frequency -\r , there exists two values of x, x,, and x 2 that 
lie on either side of some "band-center” x 0 . Of course, x 0 corresponds to d-c on the 
\-axis. 
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It can be readily shown that 


I 


x o = X 1 x 2 (78) 

From Figure 12a, it can be seen that the points x 0 , x lf and x 2 map to some band-pass 
region on the w-axis. In addition, there are other band-pass regions as expected. 

From Equations (78) and (76) the following relationship is easily derived: 

w„T w T w T 

tan —= tan-- - tan — - (79) 

Letting \ = 0, and w= w 0 in Eq. (75) results in 


tan 2 



Thus Eq. (79) reduces to 


w T wT 

tan -- - tan = a 2 

2 2 


(80) 


(81) 


IX. PRACTICAL "BAND- PASS" DIGITAL FILTER DESIGN 

As a prelude to explaining the steps in a band-pass filter design, it will be assumed 
that the following parameters are specified: 

lower "cutoff frequency = 

upper "cutoff" frequency = w 2 

bandwidth = w 2 - 

In some cases it may be desirable to specify the "center" frequency ff 0 , although it is 
more convenient to use the upper and lower frequency references, and this is the only 
case that will be considered. The design steps follow: 

1. Obtain an analog low-pass transfer function of the type desired in the band-pass 
of the digital filter from an appropriate approximation technique. The most con- 
venient choice of "cutoff frequency in the low- pass prototype is \= 1. 

2. Determine the constants a 2 and A from the relationships: 


w, T wT 

tan—— tan 

2 2 


(82) 


A =- 


\r 

w 2 T w.T 

tan-— — tan 

2 2 


(83) 


where \r is usually the reference low-pass "cutoff frequency. 
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3. The required transformation is 


\ 


P = A 


I" 2 -1 

Lz + 1 + (z-l)j 


(84) 


This expression is substituted into G<p) thus yielding H(z). 

4. The quantity H(z) is synthesized according to the techniques of Sec. m. 

5. K desired, the frequency response of the filter may be plotted with the aid of 
Figure 13 and Figure 8 of Sec. VI. 



X. BAND PASS DESIGN EXAMPLE 

As an illustration of the design procedure, a band-pass filter will be designed using 
the low-pass prototype of Example 2 in Section VU. The low-pass prototype was a 
Chebyshev filter with 1 db ripple in the pass band and a normalized cutoff frequency of 
\ = 1. Suppose the following band-pass specifications are given: 

f s = 2 kc 

f j = 300 cps = 0.3 kc 
f 2 = 400 cps = 0.4 kc 
BW = 400 - 300 = 100 cps = 0.1 kc 
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The low-pass prototype is 


G(p) = i (85) 

1 + 2.6943p + 5.2748p 2 + 3.4568p 3 + 3.6280p 4 

The constants a 2 (along with a) and A are determined from Equations (82) and (83) 

a 2 = tan .15^ tan J2n 
= tan 27° tan 36° 

( 86 ) 

= 0.50953 x 0.72654 
= 0.37019 

a = 0.6084 (87) 

A = — -* i (88) 

tan 36° - tan 27° 0.72654 - 0.509K3 

= 4.6081 


The required transformation is 


p . 4.6081 [ i-4 (89) 

L z + 1 (z - 1) j 

which may be inserted into Eq. (85) to yield H(z), the desired transfer function. The 
procedure is straightforward, but the computations are rather lengthy and messy for 
this case. Consequently, the actual synthesis will not be carried to completion for this 
case. As a final exercise f >r this example, the frequency response will be plotted. Fig- 
ures 8 and 13 and the low-pass tabulated data may be used in this context. Consider the 
normalized forms of Equations (76) and (77) which are repeated 


(90) 


where x = (91) 


These equations are further written as 



h—k 

(92) 

and 

ax = tan wT/2 

(93) 
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Thus the vertical axis of Figure 13 is interpreted as \/aA and the vertical scale of Fig 
ure 8 is interpreted as ax, The quantity a A is numerically equal to 2.804. Hence, 


A ___ 1 

2~804 ~ X ”'x 


and 0.6084 x = tan wT/2 

The results are tabulated below, and a plot is shown in Figure 14. 


X„ 
608 X 

a 

608 X,, 
: a /l kc 
L/l kc 


0 

.383 

.707 

.924 

1 

1.05 

1.58 

2.1 

2.63 


.891 

1 

.891 

1 

.891 

.707 

.603 

.016 

.006 

0 

0 

.137 

.252 

.33 

.357 

.375 

.564 

.75 

.94 


1 

.935 

.88 

.84 

.835 

.83 

.755 

.7 

.64 

0 


1.08 

1.14 

1.19 

1.2 

1.215 

1.32 

1.445 

1.62 

1.44 1 

.608 

.568 

.535 

.511 

.507 

.504 

.459 

.425 

.389 

0 


.657 

.693 

.723 

.73 

.739 

.802 

.876 

.973 


.35 

.33 

.313 

.3 

.299 

.298 

.275 

.255 

.233 

0 


.369 

.387 

.4 

.401 

.402 

.431 

.458 

.492 

1 


Note: Table is accurate only to s!ide~ru!e and curve accuracy 


8th Order Chebyshev Band Pass 
€ - \ db 

f s - 2kc 
fc\ =0.3Kc 
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XI. SUMMARY 

Digital Filters provide a means for altering the spectral representation of sampled 
data in the same manner as RLC filters do for continuous analog data. Furthermore, 
the versatility of the computer allows the filter characteristic to be quickly changed to 
conform to rapidly changing characteristics of telemetry data. 

In a sense all digital filters have periodic frequency response characteristics. How- 
ever, if the sampling rate is chosen sufficiently high and if accurate digit - 1 to- analog 
conversion is employed, the periodic effect of the response may be removed. 

The procedures presented in this paper depend on the bilinear transformation to 
map an analog domain response into the digital domain. The approximation problem is 
first performed in the analog domain U3ing conventional network synthesis tables, and 
then this response is mapped into the digital domain. The use of the z-transform results 
in a numerical algorithm that can be easily programmed on a computer. 
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I. INTRODUCTION 

Tracking of satellites and stars with optical devices requires that these devices be 
precisely pointed in space. A system for obtaining a high degree of pointing accuracy is 
currently being developed at the Goddard Space Flight Center by the Automatic Control 
Group of the Antenna Systems Branch. Called the "multi- mode mount," it consists of a 
two gimbal structure capable of positioning optical equipment weighing 3,500 pounds and 
having maximum dimensions of 14 feet in length and three feet in diameter. The angles 
of die gimbals’ axes correspond to X and Y position coordinates of the object being 
tracked. The specified position accuracy is 0.0001 degrees in both coordinates, and the 
system is required to move smoothly at rates from 0.091 degrees per second to 5 degrees 
per second in order to track stars as well as low altitude satellites. 

The control system for the multi-mode mount is to be a hybrid (digital -analog) type 
consisting of a "tight" inner rate loop and a relatively "loose" outer position loop. This 
same scheme is applied independently to both the X and Y axes. A functional block dia- 
gram for the X axis control system is shown in Figure 1. Digital position and rate feed- 
back signals are obtained from digital encoders located on the axes of the gimbals. These 
encoders provide position information accurate to within ±0.001 degree. Analog rate 
signals are also available from dc tachometer generators mounted on the axes. 

Though both loops operate on digital signals, the equivalent sampling rate of the inner 
loop under most operating conditions is high enough to assume that the signals are analog. 
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This is not true, however, in the ca the outer loop. Thus, a continuous signal analysis 
is applicable to the design of the inner loop, but discrete signal analysis must be applied 
to the outer loop. Quantization effects were not included in the initial analysis. Errors 
due to quantization should be small since quantization is in steps of 0.0001 degree. 

The study performed by this author relates to system input requirements m the pro- 
gram mode of operation, analysis of experimental results obtained in the rate mode of 
operation, and compensation requirements for the rate loop. In the prc^ ’am mode, the 
input is to the position loop and consists of a digital signal corresponding to the sampled 
and held value of the satellite's position T seconds in the future. Sampling is performed 
once every T seconds. The present mount position is subtracted from this signal, and 
the difference is sampled and held for T seconds. A rate signal generated on the basis of 
the sampled difference drives the rate loop at a rate required to reduce the difference to 
zero in T seconds. The rate signal is constant over the T second interval. In the rate 
mode, the position loop is open and a command rate signal is applied directly to the rate 
loop. Experimental results referred to for this mode of operation were limit cycles of 
unacceptably large amplitude which occurred when driving the system at low rates. 

The analysis of the system input requirements in the program mode showed that in 
the case of a satellite with a 100-mile orbit on a worst case (horizon) pass, the sampling 
period, T. must be less than 0.282 second to obtain errors less than 0.001 degrees. An 
alternative solution is proposed in which the input to the position loop is derived from the 
satellite's present position, velocity, and acceleration. With this input a sampling period 
of one second would be sufficient to maintain 0.001 degree accuracy. 

From an analysis of the experimental results, it was determined that the most likely 
cause of the undesired limit cycles is the stick-slip phenomena associated with bearing 
static friction. To eliminate this problem, the damping of the rate loop must be increased. 
This in turn requires a reduction of the rate loop bandwidth and a concomitant decrease 
in response time. A design was suggested in which a feedforward path in the rate loop is 
used to maintain required response time while allowing a reduction of loop bandwidth and 
an increase in loop damping. 

This report is divided into four sections exclusive of Introduction and Conclusions. 
Section II deals with pulse models of the system to which z transform analysis may be 
applied. In Section in, requirements for the reference input are discussed. An analysis 
of experimental results is given in Section IV, and some design suggestions are included 
in Section V. 


II. SIMPLIFIED MODEL FOB DISCRETE TIME ANALYSIS 

In order to arrive at a mathematically tractable pulse model for the system, some 
simplifying assumptions must be made. First the quantizing effect due to the D/A con- 
verter in the digital tach loop is neglected. Since the quantization level is quite small, 
this should not lead to serious errors. Next, the digital tach loop transfer function is 
assumed to be simply 1/s. This assumption neglects dynamics which should have a short 
transient relative to the sampling period T. On this basis the pulse model in Figure 2 
applies. 

In the figure, H and G are the Laplace transforms of the hold and digital tach loop 
respectively, T is the sampling period, t\ is the command rate, i.e., the output of the 
rate multiplier, and R(s> is the reference input which is discussed in detail in Section ID. 
The actual dynamics and design of the digital tach loop are considered in Section V. 
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Figure 2— Pulse Model for Simplified System 

Ihe basis for the z transform analysis to be applied to the pulse model may be found 
in many texts. Several are listed as References 1 and 2. The analysis is as follows: 

Let the integration of the zero order 1 hold be combined with the integration represent- 
ing the digital tach loop to form a single transfer function G t which is 

G (s) =i-G(s) =-1- (1) 

1 s 


In this way, the problem of finding HG(z) is reduced to finding the z transform of Gj (s). 
This is 


( Gj (s) } s — — 

\ 1 I <*-1 


Use of (2) allows HG (z) to be determined as 

np f .v K(z-l). Tz . KT 
HG(z) - - — 

The closed loop pulse transfer function relating & (z) and R(z) is found to be 


For KT = 1, (4) becomes 


UJz) _ HG(z) KT 

R (z) 1 + HG(z) " Z - 1 , KT 


=JL; (KT = 1) 

R(z) z 


It is necessary to consider the nature of the reference input, r (t) , to determine 
whether or not tl « closed loop pulse transfer function given by (5) is desirable. The ref- 
erence input presently used is 


I 







r (t) - 0 T (nT + T) ; nT = t SnT + T 


( 6 ) 


where f-\ (t) is the computed value of present satellite position. 

Since the reference input used is the satellite position one sample period in advance 
of present time, the following z transform relationship holds 

z (r<t)} = zZ {<? r (t>} (7) 


If (7) is combined with (5) the result is 


g(») 
*,<*> " 


( 8 ) 


where d T (z) is the z transform of 6 t (t) . 

According to (8), (5) is a desirable transfer function since with the r (t) used, it re- 
sults in the output being identically equal to the computed satellite position at least at 
the sampling instants. 

The above analysis is only approximate since the digital loop transfer function is not 
a pure integrator as assumed. It is well to investigate the situation when a better approxi- 
mation to the digital tach loop transfer function is used. Toward this end, let G(s) be 
given by 


7 


1 


Then Gj <s) is 




a, !)•]-■ 


and the z transform of G t ( S ) becomes, 1 

gl (T)z(z-/3 t )(z-/3j) 


<*!<*> =- 


(z-l) J (z-aj) 2 


where a t = e"“ lT ) J3 1 and are solutions of 


_ 2(2a ia .T + a 2-1) 

P — ^ 1 — fii a, = 0 

2a +a aj T-2+ajT 1 


(9) 


( 10 ) 


(ID 


and 

T(l-a,) 2 

gl(T) = (l-/J l )(l-/3 J ) 
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The values of and /3 2 may be found from curves on page 79 of Reference 1. 
The solution for HG (z) in the case becomes 

K gl (T)(z-^)(z-^) 

Hu(z) = 

(Z-lXz-ttj ) 2 

The closed loop pulse transfer function corresponding to (12) is 


6 Kg 1 (T)(z-/3 1 )<z-/* 2 ) 

-~( 2 ) ~ ~ 

R (z - 1) (z -a,) 2 + Kg, (T) (z -/?,) (z -/s 2 ) 


(13) 


Making use of (7) in (13) gives 

± (z) Kg,(T)z(z-^ )(z -^ ) (1 4 ) 

(z-lXz-aj) 2 + Kg, (T)(z-^,)(z-/J 2 ) 


The design parameters available in (14) are K, T, and a,. If the desired response 
cannot be obtained through selection of these parameters, then digital compensation must 
be added. This may be done either within the position loop or before it in the form of a 
prefilter. 


m. REQUIREMENTS FOR THE REFERENCE INPUT 

In order for d(t) to track 0 t (t) within a few seconds of arc, the reference input 
r(t) must meet certain requirements. It is the purpose of this section to determine what 
these requirements are. 

The system is represented Ln simplifed form in Figure 3. 

The command rate 9 C is generated from y (t) and is constant over period T. For the 
purpose of determining requirements on the reference input, assume that the loop L, is 
such that £(t) 3 h c , i.e., it follows the command rate exactly. 


m 



e 



Figur* 3-SimpIi{.»d Functional Diagram 







The command rate, and thus A, during th<~ n th sampling period is given by 


i (t) -b c (t) - y - L imJl . ? - I2II ; nT - t ^ (n + 1)T 

If the reference input is given by 

r (t) = B t (nT + T); nT ^ t * (n + 1)T 


then an equivalent representation for the system input is as shown in Figure 4. 



This representation is derived through use of the relationship 

y (t) = 0 t (nT + T) - 6 <t) = 0 t (nT+T) - 6, (n%) + (nT) - 0(t) 


Substitution of (17) into (15) gives 


0(t) = 


0 t (nT + T) - 0 t (nT) 0 t (nT) - 0 (nT) 


The first term on the right side of (18) is an approximation to 6 t (t) , while 
is an additional component of velocity to correct for the error e (nT) existing a 
ginning of the n th interval. Since the first term is only an approximation to <> r 
the ultimate error will depend on the accuracy of this approximation. In order 
onstrate this, consider an interval for which 

e (nT) = 0 t (nT) - 0(nT) = 0 


Then integration of (18) for the interval of interest gives 


0(t) = 0(nT) + 


[0, (nT + T) - £ (nT)] r 
T 


for 0 S | t where r = t - nT 
’2 


(15) 

(16) 


(17) 

(18) 

i the second 
t the be- 
(t) , then 
to dem- 

(19) 

( 20 ) 


1 
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Thus, during this interval, 

e (t) = 0 t (t) - e (t) 


where o' t (t) = 6 t (t) - & t (nT) 
It is clear from (21) that 


since by definition 


6 (nT + T) = 0 


6>' r (nT + T) = 6 t (nT + T) - 0 r (nT) 


( 21 ) 


(22) 


The maximum error occurring during the interval is at t = tj for which 

*/t = t, = 0. 


de 

dt 


Time is found from the solution of the equation 

X I . 

d^ r (t) e t (nT + T) - e t (nT) 


dt 


In general 


9; (t) = f 0 f dt; nT^t *(n + l)T 

•'nT 


Let 0 t be given by 


0 (r) = a + br 


(23). 


(24) 


(25) 


where a and b are the satellite velocity and acceleration evaluated at t = nt. Use of 
(23), (24), and (25) gives 


where A * Q t (nT + T) - 0 t (nT) 


d(? t (t) 

■ nr 3T' 


= a + br =- 


. (26), 

* , , * . I 

1 - ‘ . '! '■ : 1 < 
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Thus, the time at which the maximum error occurs is found from (26) to be 



(27) 


and since A = aT +. (bT 2 )/2, this is 



(28) 


The maximum error is 


I 

f 



(29) 




For a satellite in a 100- mile orbit, b may be as large as 0.1 degree per second 2 . 

In this case then, e nax = .0125° at t/2 = 0.5 if T is one second. Under the best cir- 
cumstances, then, this is the ultimate attainable accuracy if r(t) is chosen as in (16) 
and T = 1 second is used. Reducing T to 0.5 seconds decreases e aax to .003125°. In 
order to maintain |e aax | < .001 degree, T must be less than 0.282 seconds. 

There is a problem in trying to reduce the error by simply decreasing T. This is 
that the approximation that e = <9. will be less valid over an interval as T decreases. A 
possible way to decrease the error while maintaining T = 1 second is now considered. 

Since the error results from the fact that the command velocity differs from true 
satellite velocity over an interval, a way must be found to bring the two into closer 
agreement. One possible way for doing this is to generate the command rate from a 
function of the satellite’s present position, velocity and acceleration. The reference input 
required for accomplishing this is 


r(r) = 5 r (nT) + C‘ t (nT) + S f (nT) 


(30) 


for 0 = r = T, where r = t - nT. The command velocity would then be generated ac- 
cording to the equation 

0 f (nT)-«?(nD . 

6 c (r) = — + 0/nT) + 6 r (nT) r (31) 

An alternative method for obtaining the same result indicated above is possible even 
if velocity and acceleration signals are not available. This is to generate the command 
rate according to the equation 


h y(nT - T) . 2 ty(nT) - y(nT - T) ] r 
» C < T > f + 


(32) 


where y is defined as in Figure 3. In this case the reference input to be used is that 
given by (16). 

Either of these two approaches should give errors less than .001 degrees without re- 
quiring sampling periods less than one second. 
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IV. ANALYSIS OF EXPERIMENTAL RESULTS 

Experiments on the multi- mode mount revealed the existence of limit cycles while 
driving the system at low rates. The cause of this phenomena is most likely the bearing 
friction. A detailed exposition of limit cycles of this type may be found in Reference 3. 

A brief discussion is given below to explain their occurrence in this system. 

The analogue tach loop is shown in Figure 5a. The nature of the bearing friction, 
f B ( 6), though not known exactly, is probably as shown in Figure 5b. This characteristic 
is simplified in that hysteresis effects are neglected. Though it would be necessary to 
include hysteresis effects to obtain an accurate quantitative description for different 
drive rates, the simplified representation of bearing friction is adequate to explain the 
presence of the limit cycles. 

Consider the input to be a constant velocity, 8 i . A describing function can be derived 
for the bearing friction about f^, where 

(33) 


This describing function defines an equivalent gain for f B as a function of Ae. it is 
developed by assuming a sinusoidal variation of velocity, (hb) sin «t, and calculating 
the fundamental component in the output waveform which will vary between f dt and f B2 . 
Note that it will be a function of < 9 . as well as he . Let this describing function be 

0(0., h0) (34) 


With this definition, the effects of bearing friction can be introduced into the system block 
diagram as shown in Figure 6. 

It is important to note that the effect of the feedback of e through N is destabilizing, 
i.e., an increase in e causes a decrease in retarding torque. This is in contrast to 
feedback through D and K T , which is stabilizing. 

The transfer function relating 0 to 0. is 


a(s) 

0 t (s) 


kj K 

1 + kj KK t 

it' T + T* 

g hi 2 g m 4 

1 + k, KK t + 1 + kj KK t + 


(35) 


where k, = 1/(D -N) and r' m = (r n ) [(D/D - N)] , 

Since k t and rj are dependent on N, it is clear that this transfer function can change 
considerably with For example, from (35) it is seen that 


D - N + KK t 

r r D 

e m 


(36) 
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Thus, the natural frequency, u> , decreases for increasing n . The damping ratio may 
increase or decrease depending on relative magnitudes of other system parameters com- 
pared to N . 

The gain of the transfer function given by (35) is 


k »* K 

1 4 kj D - N + KKj 


(38) 


This gain increases for increasing N . Thus, the combination of increasing gain arid de- 
creasing natural frequency can lead to conditions reqtiired to sustain a limit cycle in the 
closed digital tach loop which has the open loop transfer function 


H2(.6S + l)» 
s J (3s + 1) 



(39) 
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One limit cycle frequency which occurred in the experiments was 5 cps when driving 
at 0.1 degree per second. It can be seen from a Bode plot of Lj (s) (Figure 7) that the 
complex conjugate poles of n V(s) could move to a position to give 180° phase shift and 
0 db gain at 5 cps. 

When driving at 0.04 degree per second, the resulting limit cycle was not sinusoidal. 
To quantitatively account for this case, hysteresis and the exact form of f B (0) would 
have to be taken into account. 


V. DESIGN SUGGESTIONS 

The limit cycles occurring at low drive rates may be substantially reduced or elim- 
nated if the response of the tach loop is made overdamped. This fact is noted in Refer- 
ences 3 and 4. Reference 4 deals with nonlinear compensation schemes intended to 
overcome this type of limit cycle. The design discussed herein does not use this approach 
since it would require extensive modification of existing equipment. Instead, the design 
suggested herein for the digital tach loop requires only linear compensating elements, 
and requires only a minimum of modification to existing equipment. 



<*0t 0.! i.o too 1000 

u 

Figure 7-Bod* Diagram for L } (s) 
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The design philosophy will be considered before discussing details. This philosophy 
is to relieve the digital tach loop of the entire burden of meeting steady state accuracies 
to sinusoidal inputs for "worst case" satellite passes. This is accomplished by using a 
feedforward path as shown in Figure 8. With this arrangement, the accuracies mentioned 
above can be obtained with much lower digital tach loop bandwidth than is presently used. 



Figure 8-Feedfcrword Compensation 


In Figure 8, d e is the command rate into the digital tach loop (shown in Figure 3), 
and it is. assumed that the closed analog tach loop has the transfer function, 



(40' 


i.e., the analog tach loop gain is high enough so that the complex conjugate poles result- 
ing from motor and electrical time constants are of high enough frequency to be neglected. 
More will be said on this point below when details of the design of the analog tach loop 
are discussed. 

The transfer function relating b to i> e is 



1 + T3 + 



*T 


fS f 1 


(41) 


This transfer function can be made identically equal to one by choosing 

K f ,K t 


(42) 


t 

t 
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To implement this feedforward loop, an additional D 'k converter is required to con- 
vert b c to an alalog signal. Note that «? c is available in digital form as (y(nT)VT (see 
Figure 3) before it is converted to pulses by the rate multiplier. 

An important consideration in this design is the saturation level of the power ampli- 
fier. It should be noted that if k t = K r = 0.715 volts/deg/sec., then an incremental in- 
crease A.-) c equal to .05 degrees per second, which corresponds to the largest expected 
value for a 100- mile orbit satellite, will produce an incremental step voltage of 35.75 
millivolts into the closed analog tach loop. This must be borne m mind when designing 
the analog loop. 

Figure 9 shows a Bode Diagram for the open loop digital tach loop transfer function. 
It includes the double break at = 40 which represents the natural frequency of the 
closed analog tach loop. This part of the transfer function was not included in Figure 8. 
Values which pertain to Figure 9 are: 


1^ = 0.715 

Kj> = 10 

K = 0.645 

C 

r = 0.477 


If the poles at «.• are neglected, the characteristic equation of the closed digital 
tach loop is given by 


S J + 4.3S + 9 


(43) 



0.01 0.1 IO 100 100.0 

w 


figure 9- Bod* Diagram lor Digital Tach Loop 
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Though not underdamped, this loop has a damping ratio of 0.715, which will produce very 
small overshoots. This damping should be adequate to avoid limit cycles at low drive 
rates. In this regard, the feedforward compensation appears to have the following ad- 
vantages. The command rate, t> c , is applied, thi nigh K f , immediately to the input of the 
analog tach loop. Even for low drive rates, this results in a signal which is adequate to 
overcome bearing friction. Thus, the mount begins to move immediately rather than 
having to wait for adequate signal level to develop in the counter as it must if feed forward 
compensation is not used. 

Attention is now directed to the design of the analog tach loop. One of the difficulties 
associated with design of this loop is that problems may arise if the load inertia, J , 
varies during operation. To pit these remarks on a quantitative basis, the closed analog 
loop transfer function is given by 


K K 


a m 


V 





j 0 fr " 


+ T e ) S + 1 


(44) 


where v is applied voltage, K a amplifier gain, K m motor torque constant (lb-ft/volt), D 
axis damping (lb- ft/deg/sec), r motor time constant, t. electrical time constant, and 
K,. tach constant (volts/ deg/ sec J. For r n = 11.8, r, = .035, K m = 6.25, D = 8.95, and 
k t = 0.715, the amplifier gain required can be calculated from 


= 1600 = 


D + K K K. 

T a m 

Dr T 


(45) 


Thus, 


1600 ( 8 . 95 ) (. 035 ) ( 11 . 8 ) - 8.95 
( 6 . 25 ) <'. 715 ) 


= 1325 


With this value of amplifier gain, the power amplifier saturates for inputs of 56/1325 = 
42.3 millivolts. Since this is greater than the incremental input of 35.75 mv occurring 
for "worst case" passes, the power amplifier should normally operate in the linear 
region without saturating. 


If r m = 11.8 is tiie largest value occurring during operation, and load inertia changes 
are not too rapid, then ^ l 40, hence the double break at ^ in Figure 9 should always 
be sufficiently far removed from the 0 db crossover frequency to be negligible. 


VI. CONCLUSIONS 


3 / iff 


A pulse model suitable for Z transform analysis of the system was developed. It 
shows that system operation is exactly as desired when secondary effects are neglected. 
The analysis of reference input requirements showed that in order to satisfy error re- 
quirements, inputs must be applied at least once every .232 seconds or, alternatively, 
rate and acceleration inputs must be used. Bearing friction was revealed as the most 
likely cause of the limit cycles observed while operating in the rate mode. A design was 


I 
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presented which should overcome these limit cycles. It was based on increasing the 
digital tach loop damping, lowering its bandwidth and incorporating a £e< 

This design requires only minimum modification of existing equipment. 
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ABSTRACT 


Gain and noise considerations make it desirable to increase the 
frequency of operation of satellites. This increase in frequency places 
more stringent requirements on the accuracy of tracking systems. 

The object of this study is to find the limitations of the 85-foot 
antenna tracking system at Rosman, North Carolina, and to suggest 
some possible modifications to the system in order to meet a specified 
accuracy of 0.01 degree. This accuracy which corresponds approxi- 
mately to a system frequency of 8,000 Me snould also u-t attained for 
low altitude satellites. 

In the absence of the exact transfer function for the system under 
consideration, a linear system was assumed. Although the model does 
not exactly represent the actual system, some insight as to the quan- 
titative values of errors was gained from its use. 

The error caused by satellite motion was calculated. It was found 
that this error is not negligible for low altitude (100 miles) satellites. 
A reduction of this error by a change in gain is not attainable without 
causing additional problems. This is due to the low resonant frequency 
of the antenna dish and the consideration of possible nonlinearities in 
the system. 

A signal adaptive technique is proposed to reduce the tracking 
error. This method introduces a computer computed function at the 
system's input. However, this requires an accurate knowledge of 
satellite velocity and position. 

The effects of wind torques and noise disturbances were also con- 
sidered. The mean square error due to wind torques and noise were 
determined. It was found that in the presence of wind and present sys- 
tem does not satisfy the accuracy requirement it should ultimately 
possess. Compensation schemes were used and a way was found to 
improve system performance of application of proportional feedback. 

As the effect of noise in the system error is dependent on the 
magnitude of the noise spectrum, a limit may be set in order to satisfy 
the required accuracy. 

It is also shown that an adaptive system could reduce the error due 
to the noise. 
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BLOCK DIAGRAM AND TRANSFER FUNCTION OF 
SERVO DRIVE AND DISH 

As the first problem encountered is the determination of tracking accuracy for the 
present system, its transfer function must be obtained. Unfortunately at the present time 
the exact system transfer function is unknown and some approximations must be made in 
order to describe the system. The structure of the dish, for instance, is quite compli- 
cated and the search for its transfer function presents, in itself, a sizeable task; however, 
it can be argued that a simplified version of the actual system function may throw much 
light upon the quantitative values of the errors in tracking. 

Figure 1 shows the block diagram of the assumed system used for this study. 



Figure 1 


A more detailed block diagram of the hydraulic system is given in Figure 2. The 
transfer functions G 2 and G 3 may be combined yielding 


K s S(aS 3 + bS 2 + eS + 1) 


Assuming the roots of the third degree polynomial in the denominator of G s to be away 
from the natural frequency w 2 of the dish, the block diagram is simplified further. Also, 
in order to introduce wind torque disturbances into the system, a block is added resulting 
in the block diagram of Figure 3. 

The simplified transfer function for the system as obtained from Figure 3 is then 


0 h K(1 + ST t ) 

E ” S 2 (1 + ST 2 ) (S 2 /a? + 24/.J, S + 1) 
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where, for the X-axis 


T, = 0.687 see 
T 2 * 0.167 sec 
W n = 13.8 rad/sec 
K = 4.84 
0 = o.l 

thus 

G(S) , 4.84(1 + 0.667S) 

S 2 (l + 0.167S) f — + 0.0145S + ll 

L (13.8) 2 J 

The open and closed loop magnitude plots are shown in Figure 4. From these plots, it is 
seen that an overshoot of forty percent is possible if the input to the system is a step 
function. Although this is quite a high value, it may not be of importance as the antenna 
in the tracking mode may not be ever subjected to such an input. 
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ERROR CAUSED BY SATELLITE MOTION 

This error is a function of satellite orbit and antenna position with respect to the 
satellite. Larger errors are present when tracking low altitude satellites than high 
altitude ones. Also, if the antenna is tracking two satellites of the same orbit height, the 
antenna error varies depending on the position of the orbit with respect to the antenna. 

For an X-Y antenna, this error would be the largest on a horizontal passing rather than 
an overhead passing as is the case for an azimuth elevation antenna. {See also the attached 
computed position, velocity, and acceleration of the antenna output.) 

Two methods may be used to evaluate the error and the antenna’s capabilities. The 
sinusoidal method which is suitable for design purposes assumes the variation of position, 
velocity and acceleration of the antenna to be sine waves superimposed to the actual com- 
puted curves. From this sine wave the corresponding position is computed and as the ac- 
ceptable error is also known, the value <? /E may be calculated. One numerically computed 
value corresponding to the 100- mile satellite is shown \ v. Also Figure 5 shows how 
the sine wave is superimposed. This calculation is for an error of 0.01 degree. 



Figure 5 

Position Curve 

T -- 160 sec 

W : - r 0.0392 

160 

111 r,< ‘ 

| = -LL 64 
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Velocity Curve 


T - 80 sec 
W = 0.078 

6 = _i JL- = 19.3° = 0.336 rad 
v 0.078 

10 4 = 19.2 x 10 2 

65.5 db 


Acceleration Curve 


T = 48 sec 

W = 2- = 0. 13 
48 

= -2££Lx 5.5° - 0.0955 rad 
(0.13) J 

£ = M2|i x 10 4 = 5.45 X I0 2 

E 1.75 

54.5 db 

With this method we plot the open loop transfer function of the system for a given 
error. The actual transfer function should have a larger gain if a lower error is re- 
quired. Figure 6 is the open loop plot of the control system. It shows that with the 
present system the error for a 100-mile satellite is larger than 0.01 degree. For a 600- 
mile satellite the error is less. Also, it is shown that for a 100- mile satellite the error 
is about two minutes. 

In Figure 7 an attempt is made to reduce the error by the change of open loop char- 
acteristic. The result is not good because of the following facts. A-60 db/dec ''lope 
should not be permitted because any nonlinearity could subject the antenna to violent os- 
cillations. Also the crossover frequency comes closer to the resonant frequency and 
consequently makes the settling time very long. The addition of zeros near the complex 
poles was tried in order to overcome this difficulty. This was not successful as the 
maximum gain for stable behavior was below the one obtained from Figure 6. Figure 8 
is an example of this attempt for 

C K(S + 1.5) (1 t 385S + 200S 2 ) 

E S(S 2 i 3S . 220) (0.07 + S) <S 2 f 20S + 400) 

It should be borne in mind that the assumed transfer function, from which the open 
loop plot was obtained, might not represent the system accurately. For this reason the 
actual error could well be more or less than the calculated values. It is hoped that ex- 
perimental field data from the present system will result in an accurate system function 
description thus eliminating the uncertainty in the calculated values. 

One possible method of tracking error correction would be the use of an adaptive 
signal according to the following: If a simple loop as shown in Figure 9a is considered, 
one may write 
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If one wants c /r = 1, one should have 


as a function of time into the system, 
this basic idea may be applied to the system. However, there could be one practical 
drawback: tne calculated error must be more accurate than the error desired. The set 
of curves attached to this report shows the position, velocity, and acceleration of 
antenna. From these curves it is possible to calculate the error. However, in order to 
obtain these curves one must know not only the speed of the target but also its position. 
If these curves are not known accurately, even with an accurately known transfer func- 
tion, the error may not be accurately calculated. If this method is to be used, then the 
correction should be done step by step. As at low frequency operation, a larger error 
can be tolerated at each run and a higher frequency and a better correction may be intro 
duced for each satellite. 


1 + - G[ G 2 


= 1 or 


This means that if one could inject the error 


EFFECT OF WIND IN THE MSE OF THE SYSTEM 

In absence of the exact transfer mnction for the system under consideration, the 
analytical model of Figure 3 was used. It was then established that the output angle and 
wind torque are related by the t. msfer function 

G (S) - — * 10 ~ 7 (1 i 0167S) 

" 0.167S 3 + S 2 + 3.22S * 4.84 


in the absence of an input signal. 

The mean square error of the system is then determined from the relation 


-2 


if" 

>•> L. 


G(s)G{- s) <t> (s) ds 


where $ WW (S) is the power density spectrum of the wind torque. 4 For our particular study 
the spectrum given in page 259 of reference 1 is used, i.e., 


♦ <$> = 




v t/i-s 2 


where 


/i,2C.V 0 /vJT 


v = frequency in rad/sec 

C. is the constant of proportionality relating the wind torque to the square 
(A the wind velocity. 
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v 0 = average wind velocity 

v* = RMS value of wind velocity fluctuation about the average velocity 

The range of frequencies of importance here is between 0.4 and 3 rad/sec. With v. = 

50 MPH, vftt) = 12.5 MPHand C, = 70 LB/FT/ (MPH) 2 ft 2 /rr = 2.44 x 10 9 (LB/FT) 2 . 
Thus 

72 1 r (2.4 X 10~ 7 ) 2 (1 + 0.167S) (1 - 0.167S) x 2.44 x 10 9 v ds 

2n ' } J- j<t (0.167S 3 4 . S 2 + 3-23S 4 4.84) (- 0.167S 3 + S 2 - 3.23S + 4.84) (S + v) (u - S) 


or 


-f _ 14.1 X l<r 5 V r" (1 + 0.167S) (1 - 0.167S) 

= 2iri l j<c D(s) D(- S) ds 


where 


D(S) = (0.167S 3 + S 2 + 3.23S + 4.84) (S + v) 


= 0.167S 4 + (0.167 j/ + 1) S 3 + (3.23 + v) S 2 + (3. 22v + 4.84) S + 4.84^ 


Using the tables in Appendix E of reference 1 and putting v = 1 rad/ sec, 

i = J_ r (V 4 0-^s >a>_pajs^ dg!!i0 ^ 72 

4 277 j X „ D(S) D(- S) 


With this value of i 4 , the MSE is 

= 14.1 X 10‘ s X 0.0172 r 24.2 x 10" 7 rad 2 


and hence the mean error 


e = 1.55 x 10 -3 rad = 0.089° 

Evidently, if our analytical model is correct, the system does not meet the specification 
of a maximum error of 0.6 minutes or 0.01°. 


SYSTEM COMPENSATION FOR MSE REDUCTION 

The improvement of tracking accuracy for the system is dictated by the feasibility 
of implementation. It would appear that the use of accelerometers is the only convenient 
and practical addition to be made to the existing system for error reduction. With these 
views in mind the effect of inclusion of an accelerometer, as shown in Figure 10 with 
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For the case of v - 1 rad/sec.and k a = 10, and using the tables of reference 1 

I s = 0.0168 


and the error is 

? = 0.0897° 

A computer program was coded to facilitate the evaluation of the mean square error for 
different values of k a . The results obtained indicated that the error increased mono- 
tonically with increasing value of the accelerometer constant K A . Thus, this approach 
fails to accomplish a reduction in the mean square error of the system. 

The effect of rate feedback in the inner loop of the system of Figure 10 (H = (Kg/K) s) 
was investigated next. The transfer function was found to be 

„ (g) ± 2.4 x 10~ 7 (1 + 0.167S) 

(0.167 + 12.26 x 10 _s Kg) S 3 + (1 + 7.33 x 10" 4 Kg) S 2 + 3.23S + 4.84 


and 


°* (S) . 2.4 x 10~ 7 (1 4 0.167S) 

S + v (0.167 ^ 12.26 x 10* 5 K r ) S 4 + t(0.167 + 12.26 x 10‘ s Kg) v + 1 + 7.33 x 10“ 4 J S 3 

1 

+ [(1 + 7.33 x 10‘ 4 ) u + 3.23] S 2 + (3.23v + 4.84) S + 4.84 v 

Substituting in (1) and using the aforementioned tables with v - 1 rad/sec, Kg = 10 


I 4 = 0.0173 


and 


e = 0.0892° 


Results obtained from a computer program for different values of Kg showed no improve- 
ment over the previous considered case of accelerometer feedback. 

Finally a signal proportional to the output signal was fed back in the inner loop 
(H = Kp/K). The transfer function for this case is 

n (S) 2.4 x 10~ 7 (1 +0.167S) 

0.167S 3 + (1 + 12.26 x 10” s Kg) S 2 + (3.23 + 7.33 x 10" 4 1^,) S + 4.84 
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°p^ _ 2.4 x 10~ 7 (1 + 0.167S) 

St 1 ' 0.167S 4 + (1 + 0.167v + 12.26 x 10" s Kp) S 3 + [ (1 + 12.26 x 10‘ 5 ) v + 3.23 + 7.33 x 10* 4 ] S 2 

1 

+ [(3.23 + 7.33 x l<r 4 ) v + 4.84] S + 4.84 v 


Thus, using equation 1 and the tables for K p = 10 6 and v = 1 

I 4 = 0.1395 x 10' 3 


and 


€ - 0.00802° = 0.49 


minutes which meets the specifications. 

it was found that for this case the error decreased monotonically with increasing Kp . 

The implementation of this last scheme for improving the accuracy of the system, 
although not impossible, may prove difficult. The difficulty lies on the fact that an ac- 
celerometer and two integrators are required and the bias of the accelerometer would be 
dependent on the antenna position. However, as the existing system incorporates encoders 
with a resolution of one mil at the output axes, it should be possible to program the ac- 
celerometer bias based on the encoders output. 


RMS TRACKING ERROR DUE TO NOISE 

The effect of random noise on the servo system can be evaluated as an rms servo 
error. Assuming additive gaussian white noise, the rms servo error may be calculated 
using the same scheme utilized for wind torque effects. 

The mean square error is determined from the relation 

f G„<S>G n (-S)* n (S>d S 

where 

= mean square error 

4 > = noise power spectrum at the input of the servo system, equal to a 
constant for white noise 

G n (S> = servo system closed loop function with respect to noise 
The block diagram of the system with the input noise is shown in Figure 11. 
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Figure 11 


Then 

\ iJilfV 

o K io s 2 q + ST 2 ) K„g + STj) 

N i + K 2 a + ST i > " K 10 [S 2 (1 + ST,) + K,g + ST,)] 

S 2 (l + ST,) 

and 

-p f’“ K 2 ,g + ST,)q-ST,)d S 

" K 2 0 2«j JLj« [S 2 (l + ST a ) + K 2 (l + ST,)] lS 2 (l - ST,) + K,(l - ST,)) 


The total servo gain, k,, = 4.84 and k 20 depends on the beamwidth which in turn is de- 
pendent on the frequency of operation. The expression for K, 0 is 


K 10 = 


0.7 

a 


where 

a = half beamwidth 

The mean square error equation may be written as 

I T-fhlt 4 JL f j “ C(S)C(-S)ds 
" Vio) n 2 "i I,. D(S)D(-S) 

where 

C(S) = 1 + ST, 

DfS) = T,S 3 + S 2 + K,, T, S t K,, 
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For Tj = 0.667, t 2 = 0.167, K n = 4.84 and using Reference 1 

I (0.667)* (4.84) (0.167) +0.167 

3 2(4.84) (0,167) (-4.84 (0.167) + 3.23] 


and hence 


I 3 =0.135 


So 


(4.84) 2 (0.135) 


The following table shows the error due to noise for a fixed beamwidth. This is done 
for various radio frequencies and different values of noise spectrum as the actual values 
of the latter were not known for the system considered. 


Table 1 



5.2 

0.23 

0.087 

5.2 

5.2 


13.4 

300 

800 

13.4 

13.4 


Assumed 


10*2 

i0"2 

10" 1 


10' 3 
5 x 10‘ 3 



1.76 x 10' 4 
3.52 x 10'* 
4.95 x 10' 7 
1.76 x 10" s 
8.80 x 10‘ s 


RMS 

error 

degree 


0.75 

3.4 x IQ' 3 
0.04 
0.024 
0.54 


In order to see die effect of beamwidth on die noise, we rewrite l 3 as 


(T* K n f 1) T 2 


3 K„T,[-K„T, 


* * *11 


Also the open loop transfer function may be approximated by 


G(S) = 


0.1 a;* 
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From the above two relations one finds 


e 2 ~ 0.071 <t> u> a 2 

n c 


From the table it may be seen that for some values of noise spectrum the error would 
not be acceptable. For instance a noise spectrum of 10' 2 and a frequency of operation of 
136 Mi, result in an unacceptable error. However, the above relation shows that if the 
noise spectrum occurs for a high altitude satellite, the error could be corrected by beam- 
width reduction. For instance, if the signal to noise ratio is such that the 10 _a spectrum 
applies for a 600- mile satellite, the gain needed is 25 db less than for a 10C- mile satel- 
lite. This gain reduction reduces W c by about a factor of four and the error by about a 
factor of two. From this one is led to the conclusion that an adaptive system is desirable. 


ANALYSIS OF DATA 

Field tests were performed in an attempt to establish an accurate transfer function 
for the antenna dish structure and servo drive system. The data was obtained from ac- 
celerometers placed in convenient parts of the dish, the feedbox, bull gear, tachometer, 
and other locations. 

If the accelerometer transfer function is assumed to be ks 2 , the following methods 
should help in the determination of the system transfer function from the experimental 
measurements. 

A. Linear System Assumption 

1. If the response of the open loop system to a step function is known, an ap- 
proximate frequency response for the open loop transfer function can be 
found. 

2. The cross correlation of system input and output for a white noise input 
yields the system's impulse response from which the transfer function may 
be determined. 

3. The measurement of the impulse response itself should determine the trans- 
fer function. 

4. From variable sine wave inputs the transfer function may be established. 

B. Nonlinear System Assumption 

The first step in this method is tee funding of any nonlinearities that would hinder 
the operation of the antenna system. Proust's linear method may be used in de- 
tecting such nonlinearities. In a linear system of infinite bandwidth, if a square 
wave of frequency w is applied at the input, tee Fourier series of the input and 
output are given respectively by: 

4K r I 

Input = — } sin cot + -i- sin 3<ot + 

w L 3 

Output = ~~ (oi) sin [«t + <t>(oS)} 


— sin 5o>t + 
5 


+ JL A(oj) sin [3o>t + 0(3<*>)J + 
3 


•} 
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If the linear system has a limited bandwidth, the output would not be as shown above. 
However, as the amplitude of higher harmonics are small, if the frequency w is chosen 
in such a way that only frequencies above the cutoff are not properly represented, it is 
still possible to find the useful part of the frequency response. Proust suggests that w 
should be chosen such that the twelve db cutoff frequency of the system corresponds to the 
eleventh harmonic in the Fourier expansion. If the system is highly nonlinear, the ampli- 
tude and phase of the harmonics would not agree with the Fourier series. 

Some types of nonlinearities may also be determined by additional tests. For instance 
static friction, which could be a source of instability in a system, may be detected by using 
a ramp input. 

Field tests performed have provided open loop data for square wave, impulse function, 
white noise, sine wave inputs as well as some closed loop responses. A resonant fre- 
quency of two to two and a half cycles appears to be the predominant frequency of the 
system. A damping ratio of one tenth was also observed However, some unexpected 
output variations were observed from the accelerometers’ recordings. This could not 
be explained readily. For instance, the approximate amplitude readings on the accelerom- 
eters located at the feedbox, X-bull gear and Y-bearing for a sine wave input in the Y- 
direction are given in the following table: 


Table 2 
Amplitudes 


•9 


\ 

t 


f 

cps 

X-Bull Gear 

Y- Bearing 

Feedbox 

i 

} 

i 

1 

i 

0.01 

14.5 

7.0 

3.5 

{ 

i 

0.02 

7.0 

4.0 

1.5 

i 

0.03 

4.5 

2.5 

1.2 

i 

! 

0.05 

3.0 

2.0 

- 

■t 

0.08 

2.5 

1.5 

1.2 

i 

0.10 

2.0 

1.0 

1.4 

4 

0.15 

1.5 

1.2 

- 

1 

0.25 

2.0 

1.5 

— 

i 

0.34 

3.0 

2.0 

— 

} 

0.51 

3.5 

— 

— 

\ 

0.78 

4.5 

2.5 

2.0 

* 

i 

1.00 

5.0 

3.0 

3.0 

1 


From Table H it is observed that around 0.1 to 0.15 there is a dip in amplitude in 
all three accelerometers. In the X-bull gear there was also a shift of 180 degrees in the 
phase angle at this frequency. At first it was assumed that the structure had a low res- 
onant frequency which was the cause of this amplitude variation. Figures 12a, 12b, and 
12c show the block diagram of the accelerometers and their responses. Under the above 
assumption, the change in phase angle should be as shown in Figure 12c. Nevertheless, 
the accelerometers did not exhibit this change in phase angle. 


The X-bull gear showed a sudden change of 180 degrees while the Y-bearing and 
feedbox do not exhibit any phase variation in the range of frequencies between 0.1 and 
0.15 cps. Thus the original assumption of low resonant frequency of the structure was 
not correct. 
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A second look at this situation did show that gravitational effects on the accelerometers 
could account for this anomaly. This may be proven mathematically. Figure 13 shows the 
two accelerometers located at the bull gear and the feedbox. The output voltage of an ac- 
celerometer, neglecting gravitational effects, is proportional to its ac-aleration. For 
instance, with an angular velocity o' * a cos wt, the angular acceleration is & " = - Aa 
sin wt and the linear acceleration is a = _ r Aw sin wt. However, in the presence of a 
gravitational force, and referring io Figure 12b, a force proportional to sin o is present. 
As ■ = A/ • sin wt, the force will be proportional to sin (A/w sin wt). The output com- 
ponent of the accelerometer due to such a force obviously depends on the location of the 
accelerometer with respect to the rotating axis. For a sine wave input velocity the total 
output voltage of the accelerometers at the feedbox and bull gear could then be written as 


e 


t 



sin cot 


) 


x ho 
S 


sin cot 


l 
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In order to justify the measured results, the coefficients uf the corresponding Bessel 
series were computed with the aid of a computer and for the case of the bull gear Table 3 
summarizes the results obtained. 


Table 3 


A, rps 

f 

CPS 

CO, rps 

0= A 

CO 

J(*> 

rAff 

g 

e 2 /B G 

0.035 

0.01 

0.062 

0.565 

0.27 

0.016 

-0.254 

0.035 

0.10 

0.62 

0.0565 

0.027 

0.16 

0.133 

0.035 

1.0 

6.2 

0.00565 

0.0028 

1.6 

1.6 



k = 0 


(- D k 

k! (k + n) 



This explains the reason for the change in amplitude and phase angle observed. 

The noise present in most recordings of field data made impossible the determination 
of the exact transfer functions. In particular a 10 cps frequency was very noticeable. 

As this is close to the hydraulic’s resonant frequency, it makes the presence of the latter. 

An attempt was made to analyze the data further by means of a spectrum analyzer 
with very little positive results. Probably the presence of noise of considerable ampli- 
tudes plus calibration problems and nonlinearities account for this fact. 

In conclusion, a 2.5 cps frequency appears to be the main resonant frequency of the 
system. The associated damping ratio is 0.1. It is also observed that a 12.5 cps fre- 
quency is present at some points in the feedbox and dish. 

Correlation techniques were also employed as an attempt to further analyze the field 
data. A frequency of about 1.5 cps is apparent in some of the correlograms obtained. No 
explanation has been found for this and it is suggested that further work be done in both the 
spectrum analysis and correlation ends of the data evaluation. 
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The spin- hamiltonian formalism, or more precisely, the effective spin- hamiltonian 
formalism, is a very useful one for calculating magnetic properties of paramagnetic ions 
imbedded in a crystal lattice. These properties are very significant in maser work. For 
instance, a typical calculation, which we shall later make, is that of the energy levels 
versus an applied d.c. magnetic field. Once these calculations are made one can then in- 
vestigate such things as pumping schemes, idler and signal frequencies, etc. One can 
even decide if a particular material is suitable for a maser on the basis of spin- hamiltonian 
calculations. 

The fundamental idea behind the formalism is to describe the interaction of a para- 
magnetic ion with its crystalline surroundings and an applied magnetic field by a function 
of the total spin resulting from the unpaired electrons within the ion. 

The theoretical bases for prescribing the form of the spin function and the calcula- 
tions resultl • therefrom are, of course, quantum mechanical. The coefficients appear- 
ing in toe spirt function (hereafter called toe spin- hamiltonian) ,' though calculable in theory, 
are usually determined empirically. Indeed, this forced union of experiment and theory 
is what makes the formalism so powerful. We attempt to describe the theoretical be- 
havior of a paramagnetic ion using a few emepirical constants. In this report we shall use 
two specific examples as a means of surveying the theory of the spin- hamiltonian. The 
first example is that of F +++ imbedded in Ti0 2 (rutile); the second is Ho ++ in a CaF 2 
lattice. 


I. Fe +++ : Ti0 2 [Bleaney and Stevens, 1953; Bowers aud Owen, 1955] 

The 5 unpaired electrons in Fe 3+ are each in the 3d state. This means that the ion 
itself is in an orbital s state, i.e., L = 0. To prove this we use Hund’s rule together with 
the Pauli exclusion principle [Pake, 1962] . Hund’s rules are: 

1. Assign maximum s (spin) consistent with the Pauli principle. 

2. Assign maximum L (orbital angular momentum) consistent with toe s. L is de- 
fined to be the maximum value of toe sum of the z-components of orbital angular 
momentum for the group of electrons. 

Thus, each electron has the same energy quantum number, 3, the same orbital angular 
momentum quantum number, 2 (corresponding to the d state), and, if we are to assign 
maximum spin to the electron group, the same spin quantum number, 1/2. If there is to 
be no violation of the Pauli principle, therefore, each electron must have a different 
quantum number, m , corresponding to the z-component of orbital angular momentum. Be- 
cause 1 = 2 for a d-state, we have m = 2, 1, 0, -1, -2. Thus electron #1 has m = 2, #2 
has m = 1, etc., to #5 having m = -2. The total M = m, + m 2 + m 3 + m 4 + m 5 =0. But 
since any arrangement of the five electrons among the five m-states always yields M = 0, 
we conclude that L = 0 (recall that L = u ). Thus, Fe 3+ is in an S-state (L = 0) with 
a spin equal to 5/2. 
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The spin-hamiltonian for the Fe 3+ : T10 2 complex is [Carter and Ckaya, I960] : 


H = g/3H-S + D(S 2 - 35/12) + £(S 2 - S 2 ) 


+ (a/6) (S* + S* + S« - 707/16) 


(1-D 


+ (7/36) F(S* - (95/14) S 2 + 81/16). 


The nominal values of the derived constants are: 


g = 2.0, D = 20.35 gc, E = 2.21 gc 
a = 1.1 gc, F - - 0.5 gc 

p is the Bohr magneton divided by h, Plank’s constant, and is equal to e /4-mc - 0.0014 gc 
per gauss (e = electronic charge, m = electronic mass, and c = velocity of light). H is 
the applied <tc. magnetic field, measured in gauss. Tho spin matrices S x , S and S x are 
dimensionless. Thus the hamiltonian gives us frequency directly, rather man energy. 

Of course, energy = hx frequency. 

The first term in (1-1) gives us the interaction of the total spin with the external 
magnetic field. The remaining terms represent the interaction of the ion with its crystal- 
line surrounding. To understand this, consider that the crystalline field is derived from 
an electrostatic potential which satisfies Laplace’s equation 


o. d-2) 

3x 2 3y 2 3z 2 


The solution of this equation may be expanded in a series of spherical harmonics, 

V -'Z ] d- 3 ) 

i, m 

If we consider the 1 = 2 term and use the fact that x = r sin 6 cos <p, y = r sin 6 
sin <p , z = r cos 6 , we have for the five 1 = 2 terms (for 1 = 2 we have m * 2, 1, 0, -1, 
- 2 ) 

xy, xz, yz, 3z 2 - r 2 , x 2 - y 2 . (1-4) 

If we replace the space coordinates x, y, z by the spin operators S x , s y , and S„ respec- 
tively, we get the additional terms of the spin-hamiltonian. The justification for this 
substitution lies in the Wigner-Eckart theorem of quantum mechanics [Merzbacher, 1961; 
Slichter, 1963 J. 

The x-, y-, and z-axes appearing in the spin-hamiltonian are referred to the crystal 
axes by 
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Hamiltonian axis 

X 

y 

z 


Direction in crystal 
titoj 
[ 001 ] 

[ 110 ] 


The Pauli spin- matrices, S x , S y and S x are the 6x6 matrices 



r° 

5 1/2 

0 


0 

0 

0 



5 1/2 0 

gt/2 


0 

0 

0 


1 
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gt/2 

0 


3 

0 

0 


s 










* 2 

0 

0 

3 


0 

$1/2 

0 



0 

0 

0 

gl/2 

0 

5 1/2 




0 

0 


0 

S l/2 

0 



r° 

-i 

•5 1/2 

0 


0 


0 

0 


i- 5 l/2 

0 

-i-8 1/2 

0 


0 

0 


0 

i* 

gl/2 

0 


-i*3 


0 

0 

1 

0 


0 

i* 3 


0 

-i8" 2 

0 

2 

0 


0 

0 


i*8 1/2 

0 

-i S 1 / 2 


0 


0 

0 


0 

iS 1/2 

0 


• 

~5 0 

0 

0 

0 

0~ 





0 3 

0 

0 

0 

0 




l 

0 0 

1 

0 

0 

0 



S 









X 

= T 

0 0 

0 

-1 

0 

0 





0 0 

0 

0 

-3 

0 
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0 

0 

0 

-5 




(1-5) 


We calculate, using ordinary matrix multiplication, the following matrices which 
appear in (1-1): 




“ 5 

0 

40 i/2 

0 

0 

0 — 

1 


0 

13 

0 

72 i/2 

0 

0 

\ 

S 2 = I 

40 i/2 

0 

17 

0 

72 1/2 

0 

• 

* 4 

0 

72 l/2 

0 

17 

0 

40 i/2 



0 

0 

72 1 '* 

0 

13 

0 



_ 0 

0 

0 

401/2 

0 

5 _ 
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s 2 =. 

y 



5 

0 

“( 40) 1/2 

0 

0 

0 


0 

13 

0 

-( 72) 1/2 

0 

0 

1 

-( 40)* /2 

0 

17 

0 

-( 72) 1/2 

0 

4 

0 

. ( 72 ) i/2 

0 

17 

0 

-( 40)* /2 


0 

0 

-( 72) I/2 

0 

13 

0 


0 

0 

0 

.( 40 ) 1/2 

0 

5 


S 2 =:- 


25 

0 

0 

0 

0 

0 


0 

9 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


0 

0 

0 

1 

0 

0 
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0 

0 

0 

9 

0 


0 

0 

0 

0 

0 

25 _ 


s 4 =- 


16 


and 


S 4 — 

y 16 


~ 65 

0 


44* 10 1/2 

0 24-5 1/2 0 

i 


0 

241 


0 

180-2 1/2 

0 

24*5 ,/2 


44-10* /2 

0 


401 

0 180-2 ,/2 0 



0 

180-2 1/2 

0 

401 

0 

44-10 1/2 


24«5 1/2 

0 

180-2 1/2 

0 241 

0 



0 

24*5 1/2 

0 

44-10 1/2 

0 

65 

— 


~ 65 

0 


-44-1 0 1/2 0 


24*5 1/2 


• 1 

0 

241 

0 

-180*2* /2 


0 

24-5 l/2 

-4440 1/2 

0 


401 

0 

- 

180-2 l/2 


0 

0 

-180*2 1/2 

0 

401 


0 

-44-10 1/2 

24*S 1/2 

0 


-180-2‘ /2 0 


241 


0 

_ 0 

24-5 1/2 

0 

-4Vl0 1/? 

0 

65 J 



“625 

0 0 

r - 

o 

O 
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0 

81 0 

0 0 0 






1 

0 

0 1 

0 0 0 





S 4 

A 








* 

= HT 

0 

0 0 

10 0 







0 

0 0 

0 81 0 







_0 

0 0 

0 0 625_ 






( 1 - 8 ) 


The constants 35/12, 707/16, and 81/16 are really the matrices (35/12)0, (707/16) U, 
and (81/18) U, respectively, where u is the 6 x 6 unit matrix, having l’s along the main 
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diagonal and zeros elsewhere. Using the results obtained so far we have for the spin- 
hamiltonian, H, the 6x6 matrix 



In writing down this final form for H we have_assumed that H = H a t , (a, being the unit 
vector in the z-direction), which means that H-s * HS x . 


Derivation of Energy Levels Using the Spin-Hamiltonian 

According to the general rules of quantum mechanics, die state of a system is defined 
in terms of the solution of Schrodinger’s equation: 

(1-8) 


If H does not depend on time then it is convenient to assume that 

'P = u exp (1-9) 

where u is independent of time. Upon substituting (1-9) into (1-8) it follows that u must 
satisfy the time-independent wave equation: 

Hu = Eu . (1-10) 


1 


It is apparent from (1-9) that E is an energy parameter because E/ff must be angular 
frequency. Thus, (1-10) is an eigen- value, or characteristic- value, equation in which u 
and E are the eigen-vector (or eigenstate) and energy respectively, of die system rep- 
resented by the Hamiltonian H. 

Because H is a 6 x 6 matrix, the state vector, u , must be a 6-row, 1-column vector: 



I 


Operating on (1-11) with (1-7) gives us, for (1-10): 

(0.007H + 68.15) u, + 6.98 u 3 + 1.23 u $ = Eu t 
(0.0042H - 14.78) u 2 + 9.37 u, + 1.23 = Eu 2 

6.98 u, + (0.0014H - 53.57) u 3 + 9.37 u 3 = Eu 3 

( 1 - 12 ) 

9.37 u 2 + (-0.0014H - 53.57) u 4 + 6.98 = Eu 4 

1.23 u } + 9.37 u 3 ■+ (-0.0042H - 14.78) u 5 _ Eu $ 

1.23 u 2 + 6.98 u 4 + (-0.007H + 68.15) u s = Eu 6 

or, after transposing the right-hand side: 

(0.007H + 68.15 - E) Uj + 6.98 Uj + 1.23 u s = 0 
(0.0042H - 14.78 - E) u 2 + 9.37 u 4 + 1.23 u 6 = 0 
6.98 u, + (0.0014H - 53.57 - E) u 3 + 9.37 u s = 0 

(1-13) 

9.37 u 2 + (-0.0014H - 53.57 - E) u 4 + 6.98 u { =« 

1.23 ti, + 9.37 u, + ( - 0.0042H - 14.78 - E) u = 0 

1.23 u 2 + 6.98 u 4 + (-0.007H + 68.15 - E) u # = 0 

This is a homogeneous system of linear equations to be solved for u p u 2 , u 3 , u 4 , 
u 5 , u 6 . In order for at least one of these u's to be non-zero, it is necessary that the 
determinant of the coefficient matrix vanish. We require, therefore, that: 



(0.007R+ 
68.1 SHE) 

0 

6.98 

0 

1.23 

0 


0 

(0.0042H- 

14.78-E) 

0 

9.37 

0 

1.23 

det 

6.98 

0 

(0.0014H- 

53.57-E) 

0 

9.37 

0 


0 

1.23 

9.37 

0 

0 

9.37 

(-0.0014H- 

53.57-E) 

0 

0 

(-0.0042H- 

14.78-E) 

6.98 

0 


0 

1.23 

0 

6.98 

0 

('-0.007H+ 

68.15-F.) 


For each value of the magnetic field, H, there correspond six values of E satisfying 
(1-14). A program for computing these eigen- values was written by Mr. Henry Miller 
for the IBM-7094. This program gives not only the eigen- values, E, but also the cor- 
responding eigen-vectors, u . The numerical results are shown in Fig. 1. 

It is noted that the zero- field energies occur in pairs. Each pair is called a Kramers' 
doublet to signify that it may be interpreted by Kramers' theorem which states that in a 
system containing an odd number of electrons, an electric field (in this case the crystal- 
line electric field) cannot completely remove all degeneracies (when two or more energy 
states coincide, the corresponding energy level is said to be degenerate). 
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We may calculate the energy levels for different orientations of H by modifying the 

H-S term appearing in (1-1). Figure 2 illustrates the spherical coordinates to be used. 

In terms of the spherical coordinates we have: 

H = H sin 9 cos 9 + H sin 8 sin <p a y + H cos 8 a t , (1-15) 

where 

H = I H j . 


Thus 


H*S = (H sin 9 cos <p) S x + (H sin 0 sin 9 ) S y + H cos 6 S t , (1-16) 


where S x , S y , and S f are the spin matrices of (1-5). The effect of (1-16) is to introduce 
new terms in the spin-hamiltonian matrix, H. This modifies (1-7) and (1-14) showing that 
the energy levels vary with e and 9 and, hence, with the orientation of H . 


II. Ho ++ : CaFj [Sabisky and Lewis, 1963] . 

Holmium is a type 4f rare earth, which means that the divalent Holmium ion has its 
unpaired electrons in the 4f shell where they are effectively screened from their crystalline 
surroundings by electrons in the outer shells. Therefore, as a reasonable approximation 
to the effective spin-hamiltonian we may discard any terms that represent the crystalline 
field. We must, however, include die spin-spin interaction between the unpaired electrons 
and the nucleus because these electrons are relatively close to the nucleus. Thus, we use 
the following spin-hamiltonian 


h = g/3ff •§ + aT’S, 


(2-1) 
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where g = 5.91, fi - 0.0014 gc/gauss, A = 3.924 gc, S is the electron spin operator, with 
effective spin 1/2, and I is the nuclear spin operator with spin 7/2. 

Because we are dealing with a system of two particles (electron plus nucleus) we can-* 
not simply form matrix products in order to evaluate H . In (2-1) the term T • s must be 
interpreted as the direct-product of the operators I and s. This concept will be explained 
in more detail. 


The effective electron spin of 1/2 means that there are 2 x 1/2 + 1 = 2 possible pro- 
jections of the electronic spin vector upon (say) the Z-axis: spin "up" and spin "down”. 
The nuclear spin of 7/2 implies that there are 2 x 7/2 + 1 = 8 possible orientations of the 
nuclear spin moment. Hence, there are 16 possible states jointly involving the nucleus and 
electron, eight corresponding to electronic spin up and eight more for spin down. 


If 


a, = [1 0] 


a 2 =[0 1], 


(2-2) 


n up 

and 

spin down, 

respectively, and 


4 = 

[10 

000000], 

ft s = 

[00 00 10 

0 0] 


[0 1 

000000], 


[000001 

0 0] 

P 3 = 

[0 0 

1 0 0 0 0 0], 

07 = 

[00 00 00 

1 0 ] 


[0 0 

0 1 0000], 

^8 = 

[0000 00 

01 ] 


(2-3) 


stand for the eight nuclear spin states, then the direct-products of the electronic and 
nuclear spin states are represented by the 1 x 16 matrices (or row- rectors): 

li> = [i o o] 

| 2 > = [0 1 0 0 ] 

|3> = 100 1 0 0J 

1 4> = tO 0 0 1 0 ♦ • • • 0] 


(2-4) 


etc. 


|1S> = {000 •*• • • 0 1 0] 

|l6> = [0 0 0 0 1]* 

The first eight direct-product states correspond to electronic spin up and the last 
eight to spin down. 

If we arbitrarily take the magnetic field, It, to lie along the z-axis, then = 

k ft HS t , where H is the magnitude of if. We must treat this operator as the direct product 
of s t and u, where 
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(2-6) 


is the 8 x 8 unit (or identity) matrix. 

The direct product of these two matrices is given by the 16 x 16 matrix 


S x U=4| 
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(2-7) 


Next, we must calculate the direct- product, I x §, of T and §. By definition Tx s = 
I xs x + i y x s y + l t x where, for the system being studied: 


s ‘4[! i] vi[0 S.O.J. y 


(2-8) 
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I =- 
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The required direct-products are given by the 16 x 16 matrices: 
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(2-9) 


Upon substituting (2-7) and (2-9) into (2-1), we get (2-10) as the final form of the spin- 
hamiltonian matrix. All entries not shown are zero, and H is measured in kilogauss. 

From this point on die analysis proceeds as in the treatment of iron-doped rutile. 

We mast calculate the eigen-values of (2-10) as a function of H. The results of the IBM- 
7094 computations are shown in Figure 3. 

A comparison of Figures 1 and 3 shows that Ho++ has a much more uniform varia- 
tion of energy (and hence resonant frequency) with H than does Fe ++ . This follows, as 
has been mentioned before, because the unpaired electrons in Hot* are screened from 
die crystalline field, whereas those of Fe 3+ are not. Hence, Ho 2 * behaves, for H greater 
than 3 kilogauss, as a free spin in an external magnetic field. We would expect this same 
qualitative behavior for many of the 4f rare earths, no matter what the crystal lattice is. 
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A PRELIMINARY STUDY OF HOLOGRAPHY 
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I. SUMMARY 

This is a preliminary study of holograms for information processing, attention being 
directed toward their formation, reconstruction of their images, and the images resolu- 
tion. The formation of holograms is studied with respect to information preservation for 
two- and three-dimensional objects and resolution of the interference patterns controlled 
by frequency dispersion, beam collimation, and relative motions of the laser light source 
and objects to be imaged. Here, the emphasis has been placed on contrast and energy 
ratios; so that a suitable quantization measure can be determined. For reconstruction of 
the images from holograms, the problem of "twinning” and its removal, multiple images 
spatially and/or frequency separated, and aperture limiting phenomena are explored. 
Finally, the image resolution upon reconstruction, manifesting the completeness of the 
hologram process, is found to depend on the previously discussed factors in such a way 
that additional study is required to obtain a meaningful measure. 


n. FORMATION AND RECONSTRUCTION OF HOLOGRAMS 

Information Preservation with Holograms 

A two-dimensional screen with ( N x , n ) elements per unit length is interposed along a 
collimated path of monochromatic radiation of wavelength \ from a source Q t to a re- 
cording medium £ of elements p e p, the observation point. Simultaneous with radiation 
from a source is a similar monochromatic, unobscured source q 2 radiating along a 
path oriented at ah angle o with respect to the normal to the medium’s plane, s • The 
information content 1 for the interposed screen £,, H - log 2 N x + log 2 N y , is preserved 
in the original mapping via Q t , and the interference pattern induced by the addition of the 
second source Q 2 suffers no loss of information content as evidenced by the equivalent 
field 2 generated by the "complementary" 3 screen, whose current distribution complements 
the original screen distribution. Thus, one sees the typical interference pattern for two 




Figure 1 -Formation of hologram 
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sources from which is subtracted the complementary field, and for black-white screens 
one easily notes there is no information loss. 

Collapsing Loss in Imaging with Holograms 


Allow a three-dimensional configuration of objects to be illuminated by a collimated 
monochromat of wavelength \. 


*<P) 
*< P) 

K P) 
-/-(P) 


= * 2 (P) 
= * 2 (P> 
= * 2 (P) 

= 2 cos 


+ ^ S (P) 


n 



+ >A,(P) 




- %(?)■ 


i(^ k ) 


e Ik»><l k 

Rp^ 


d tr. 


( 1 ) 


The radiation from the configuration has an equivalent distribution of current sources on 
the screen ij , hence the information content for the screen is preserved in its imaging 
on the recording medium. Because the distribution on the screen has a collapse of dimen- 
sion, a loss of information necessarily does occur in the formation of a hologram. 



Figure 2-Equivalent screen 
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P«P 


Supposing one could vary the observation point upon readout from the hologram, one 
could alter die "screen's” location, modify the shadowing or silouhetting that one sees 
and thus recover the information lost via collapse. 


Twinning and Loss of Information for Holograms 


In the formation of a hologram, one cannot avoid a difference of intensities for ^ (P) 
and </> a (P) , though they be coherent sources. For reference, <P t (P) and </- a <p> have 
been chosen conveniently to be of unit magnitude, hence we can refer to them as *i.o <P> 
and V-2,0 (T> , respectively; again, for convenience, we can choose </>. (P) = r <P) , 
r a complex constant, in order that we can refer back to the simplest case when r = 1. 
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Figure 3- The hologram scene 


image formed on the surface Z, for pc p, is of wave potential </> (P) , 

^(P) = (P) + _L.| I< a,) da t ( 2 ) 


where I (a t ) is toe equivalent current distribution for the screen z containing the pro- 
jected information from the geometrical array of objects placed between the source Qj 
and the screen, and the appropriate kernel function is for free-space. Were the objects 
removed, the screen’s distribution becoming uniform because of toe collimated light 
source Q t , the interference pattern at 2 would show toe difference of Q t and Q 2 ; upon 
replacement, the information on the screen would be preserved in the new interference 
pattern. One would find the wave potential <p(P) to be 

^(P) = </> 2 <P) + ^(P) - J T(^ G(p; <r t ) da. 

where I (a.) is now the "complementary current distribution, and in terms of the refer- 
ence solution «/>(P) is 

<KP) = 2 cos \- (1 _ D ^ 2<0 (P) - f IC^iGfpia^ (3) 


The new term (1 -r) 0 (P) prevents perfect cancellation in toe interference pattern, 

decreases the contrast ratio, and thus serves to destroy information stored in the hologram. 
This is a significant effect and accounts for "binary erasures” in toe readout process. 

Subject to the formation of a perfect hologram, in the sense that r = 1, one is 
able to reconstruct the image 0 and its conjugate O' by illumination of toe screen z with 
a coherent source q x . The conjugate image o', the twin, is out-of-focus with respect to 
0, unwanted, 4 and can be eliminated by various methods, turn of which are illustrated here. 
The first, due to Bragg and Rogers, 5 relies on the utilization of a complementary hologram 
Hj placed at twice the focal distance d from tot. ologram H, . It can be seen that the 
image O’ is cancelled by its opposite o* in the focal plane where the arrow appears. The 
second, due to Lohmann, 6 makes use of a perfectly conducting half-plane placed between 
toe hologram Hj and toe focal plane at o. The screen z with its interference pattern is 
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H, 

Figure 5-Lohmonn concelfation method 


like an antenna whose current distribution has been reduced to zero over half the aperture. 
The effect is to displace the radiation pattern toward the illuminated region, so that the 
Fraunhofer region contains no radiation from the dark space. The rays in both cases 
should contain no "twin" information, hence the resolution would be significantly improved. 

Multiple Image Reconstruction from Holograms 


Several images of objects irradiated by the same or different collimated light sources 
can be stored on a hologram for a multiplex communication system 7 - 8 or a navigation 
reference.® In the former application, several diffraction gratings are individually 
angularly modulated proportional to the intelligence to be transmitted or appropriately 
digitally coded, and the hologram formed contains the several channels of information; 
in the latter application, a Fresnel zone plate is optically spatially correlated, with x-ray 
emanations from stars in outer space, and the hologram contains their angular informa- 
tions. Upon reconstruction, one expects to recover the object information. 

When the objects are angular disjoint with respect to one another and when there is 
no loss of resolution in the process, the image reconstruction can retain the angular dis- 
jointness subject only to the confusion offered by the reading beam and twin, 10 - 11 or con- 
jugate, images. That this is the case is easily seen below. 

The reconstructed wavefronts from the hologram are "plane," the spreading over an 
angular sector o for a typical wavefront arising from errors in the current distribution 
on the hologram h- The intensity of scattered radiation in the forward direction j(p) is 
described as 


j(P> = _J_ 
(4w>* 



(°2> 


lWRjw, 

e 

Rpc^ 


-lkRs*r 2 

1 dcr dcr + 


1 


+ 


JL 

<4*) J 



*4i H* 



< a a> 


IkRjw. 
e * 

Rpoj 


-IkRpcr 

t dC7 > d£ v 

Rp^j 


(4) 
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where the first-term is the angularly- disjoint intensity for the real image and the second 
term, normally zero, is the correlation of the currents and their fields mutually ex- 
clusively. The out-of-focus ensemble, or set of twins, is not included in this description, 
and, in addition, to remove them entails some loss as to the number of objects and their 
sectors that can be stored. This is a consequence of the half-plane technique of twin 
removal. 


Image Reconstruction from Finite Holograms 


To determine the minimum size of a holo- 
gram for information storage purposes, one 
intuitively thinks 1 '• 13 of an aperture-controlled 
grating, the grating establishing orders of dif- 
fracted waves and the aperture their individual 
resolutions. The ruled lines on the holog; -cm 
constitute the grating, and they are interfero- 
metrically formed by the original and side il- 
luminations. Desirably, one wants a many wave 
length regular grating for image formation with 
irregularities less than a quarter-of-a- wave- 
length 1 4 to avoid ghost intensities. 15 

The ghost phenomena of Rowland 16 has been 
reported by Z. H. Heller 17 for a He-Ne laser, 

6943A, and an aperture-controlled grating of 15,000 lines/inch and width 2.5 x 10~ 2 cms. 
Comparing the resolution and ghosts for this situation against a similar one of a sodium 
source, where the slit width was 2.5 x 10" 3 cms., he found decided advantage for laser 
source; the corresponding situation for linear antennas is one of a few hundred linear 
elements contrasted to one of twenty elements. On this basis, one is led to the intuitively 
and satisfactory condition a > > \ with an upper bound determined by the extent to which 
coherent phase addition applies to the emerging wavefronts on reconstruction. 



Figure 6— Finite Pupil 


The situation of a limited aperture is easily expressed by 18 


<KP> -- 'VP) * A < p )> 


where the convolution of the hologram and aperture fields give rise to the diffracted field. 
Even in reconstruction a twin-image will appear due to the conjugate one at O'. The 
angular distribution of intensity is well-known for the uniformly-illuminated screen, and 
it is 


JW = 



e ik»-tliy*in4 <Jy, 


k :: constant 


i.e 


•> 


JW) * 


ka c tkK 

jjI/2 



(5) 


showing a typical Fresnel behavior of a (sin a / a.) distribution. Hie angular spread be- 
tween the two nulls symmetrically displaced with respect to the real image’s ray is <t> = 
2 sin* 1 (V«)» corresponding to twice the "space- sampling angle," and an excellent 
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reconstruction of the image occurs only if this angle <t> be very small. For hundreds of 
wavelengths of extent for the hologram, the reconstruction of the real image tends to 
reach a limit of resolution controlled by hologram irregularities and the sources fre- 
quency dispersion. 


m. RESOLUTION FOR HOLOGRAM IMAGES 

Interference Pattern Shift Induced by Frequency Dispersion 

What effect does the spectrum of the emitted, radiation from a laser have on the 
resolution property of a hologram? With simplifications in mind, one can start with a 
periodic one formed by two plane waves impinging on the screen at some now-zero rela- 
tive angle. Allowing one to be normal, the other inclined at an angle e , the undulations 
of the intensity of illumination follow a cosinusoidal law with well-known regularly- 
spaced maxima and minima. As the locations of peaks and valleys are sensitive to the 
wave-number of emission, its probability distribution affects the average peak and re- 
moves the null of the pattern. Although the spectrum is typically a quadratic functional, 
little error will be introduced if the nose of the spectrum, parabolic in shape, be fitted 
to a Gaussian error function and this be used to determine the ratio of average maximum 
to average minimum intensities. 

On the hologram itself the maxima shift in accordance with the relation 

fol-l ■ I , 

dXl I sin 6 I 


where m is the number of the maxima, m = 0 implying the origin y = 0, consequently 
the distribution p( y > is 


p(y> = 


l 


exp 


f <y-y 0 ) a \ 

V 2<7y* ; 


with 


wo X, 

sin & 


As the same applies to the minima, m permitted to assume half- integer values (± 1/2, 
± 3/2, , ± N/2), the expectations of the minimum and maximum intensities are 




exp 



du 



oj = value for • » 



• . 
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and 



The ratio of these intensities is then seen to be 



( 6 ) 


Inspection of the intensity ratio shows a dependence on the ’'noise” in the hologram, 
and one appreciates the significance of this ratio where one is asked for a measure 
suitable spatial resolution comparisons. The amount of information stored in the holo- 
gram also depends on this quantity. A related quantity is the Michelson contrast ratio 


I -i . 

q max min 

r t + I . 

m* „ T min 


and as in the case of Diederich and Lehmann, 1 ’ who studied roughness of lens surfaces 
on image formations, one finds 


C = 


A* - <r 2 
y_ 

A 2 + cr 3 

y 


(V 


Amplitude Quantization Induced by Dispersive Radiation 


Let a and (i be randomly distributed vari- 
ables taking into account the divergence of an 
otherwise collimated, monochromatic beam 
radiating from Q, and sympathetically from Q a . 
The wave potential «/«(P) for p e P induces cur- 
rents with random distributions on the screen 
£. In terms of the random angles. o x and d 2 , 
where e t = k, • ( r - r, ) and e 2 * it 2 • ( 7 - 1 2 ), 
one has 


^<P> = e 1 ** 





so that 


Figure 7~B«om Spreading 


J(P) = |0(P) | 2 = 4 cos 2 



Near the maxima, specifically when d t - e 2 * m<2w), » = 0, ±1, ±2, , the 

intensity J(P) is 


J(P> y 4 


[• 
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and for 2 = k,. j , 


7 


J(P)?4 j\ --^-(1 - k> 2 J 

with a statistical averse 

{j(P)/ = 4 [l -~(1 -k) 2 j . 

This result allows us to see a deterioration of the brightness of the hologram 2, when 
the radiation from Oj is shadowed by interposed objects. Though the deterioration is 
not rapid, the effect cannot be neglected. In sharp contrast, the minima are seriously 
affected, and they follow first 

J(P> ~ 4 ^ - ^LiSL n'j . m = 0, +1, ±2, ± ±N 


J(P) z [6/1 -k) -(2m + 1) v]t, 

second 

(j(P>/ = (a\) (1 - k ) 2 - {dj) 2(1 - k) (2m + 1) v + (2m + l ) 2 tt 2 
{j(P))?^ t (1 -k) 2 . 

and one finally sees imperfect cancellation with a ratio of maxima- to- minima statistically 
becoming 


P = 


(O 

<o 


p = 


l.—L-d-k) 2 

4 


cr\ (l-k> 2 


(8) 


The comparison of this result with the frequency dispersive relation reveals the 
quantal sensitivity of the respective phenomena, frequency and angular dispersions. The 
ratio k of beam angles is dependent only on the first-order relations corresponding to 
the * .ometry of the sources and hologram screen 2. Thus, the control of resolution for 
the ■ logram is centered on the required quantal sensitivities and the extent to which 
they are important in the face of fundamental resolution for the hologram itself. 
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Vibration and Contrast Ratio for Holograms 


The hologram screen I can move randomly relative to the (x, y) reference frame 
and thereby contribute to a loss of resolution. As unidirectional movement only can occur, 


* * 
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the displacement oi maxima and minima occur. 
The greater sensitivity occurs for the minima 
of intensity, and its expectation is proportional 
to the variance of motion. One looks for a 2 and 
finds it to be 


a 2 = tan 2 6 a 2 
y * 


Earlier, we saw the effect of frequency dis- 
persion and the induced shift of the hologram, 
so now we can easily see the minima’s expec- 
tation to be 



A \_ 4tan2 0 ff * 

Vain/ 


and the contrast ratio to be 


tan 2 9 cr 2 

1 i. 

A 2 


tan 2 9 cr 2 

1 + 

A 2 


A 2 - tan 2 0 a 2 


A 2 + tan 2 6 cr 2 


i n 

^Jraax/ + 


( 9 ) 


Energy Ratio, its Depen d ence on the Emulsion’s Properties 

For the emulsion, the transfer characteristic has been taken, for convenience and 
good approximation , to be 20 lir.ear with respect to intensity between two limits J Bax and 
J»in » *h e transmittance, or opacity, lying on the scale from zero to one but also limited 
to tj and t 0 , respectively. The second-order effect of graininess, accounting for spread 
of the images, has been included in P. Clark Jones’ analysis of information content of 
emulsions. It is now the object of the present analysis to use these desiderata in realiz- 
ing the energy ratio and resolution for holograms. 

We consider, fundamentally that the transmittance will be probabilistically, controlled 
over the screen i(s, y) in response to the wave- potential <//( y, z), and the density func- 
tion will be a continuous function thereof, that is 


dP [t(y, z)] = pt '•( (y, z)| 2 ] dydz for yeY, zeZ* 


What we shall look for is an increased variance cr 2 on the hologram arising from the 
random response of the screen i, and the simplest case to examine is the interference 
pattern for two plane waves. Here, the variance o -^ 2 is known to be normally dependent 
on frequency dispersion, beam spreading, and vibration phenomena for the laser source. 
These being independent phenomena, for instance, one would find 
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mo x 

2 



1 ] 


sin 9 

| j 



or tan 6\ 


a = 0, ±1, ± • • • , ±N 


the terms arising in the order given above and die variances determined at the minima of 
the interference pattern. No term exists for the "variance" of the dispersed beam situa- 
tion, the null losing its sharpness and becoming diffusively blurred. However, (me can 
still deal with the term "variance" as it does influence the qunatization for the hologram, 
tne energy ratio being 


P 



the minima determined by the variances. From foregoing analyses, we find 


p„ = 


wcr x 


A sin 6 


(1 -k)o 


+ 4 


tan 6 or ■ 


( 10 ) 


m = 0, ±1, ±2, ± ±N 

arid note an equivalent variance. 

Now, to account for the variation in transmittance. First, under the linear condi- 
tion's of P. Clark Jones, we note p is (me for those wave potentials outside die limits of 
the emulsion, is (t./t 0 ) for extreme potentials still outside the limits, and is the above 
number for those on the linear portion. Second, to reduce die energy ratio because of 
t / 1 , we observe that 
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•'ft 


t 2 dP [t(y, *)] 


for yeY, zeZ, 


<*•>= r 


t* pCt( |^(y, *) j 2 )] dt. 


But p(t) =P;t 0 ) S(t - t 0 ) + p(t 0 < t < t,) + P(tj) 8(t - t t ), so in the vicinity of a 
minimum, J mln - 0, 


r t ~ j^ 0 ** p (T nin )J dt 


<r 2 - t 2 

°t ~ c 0 * 


while for a maximum, J % 4, 


= <l-t,) 2 . 
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Note these are fluctuations with respect to the nulls and peaks, respectively. Clearly, 
the variance is intensity sensitive and the same energy ratio t t /t 0 is still realized under 
extreme conditions. However, in the linear portion, we find 


and 


< t2 


\ 

/ 



t 2 p [t (T)] dt 



(t -t> 2 




At 




dt 



1-P(t 0 )-P(t 1 ) 1 

At J 


dt 




' 1 -P(tp) -P(t,) ~ 
At 


(At) 3 

12 


when a uniform distribution is assumed. This variance implies that the energy ratio is 




A sin 9 


(1 - k )a e 


+ 4 


4 

tan flcr^ P ^ | "l ~P(t 0 ) -P(t t ) 

A I L 12 

ffl = 0. ±1, ±2, ±, 


At 2 J 

, ± N 


(ID 


when T min > t Q > 0. Usually At << t , the graininess, hence the simpler relation for the 
energy ratio is the one previously derived. 

As a result of the study presented above, we see that die resolution was not obtained, 
the variance not being determined readily, although the energy ratio fell on it quite 
naturally. Additional study should reveal the resolution for the hologram. 
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CORRELATION BETWEEN JITTER AND INFORMATION RECOVERY 
IN PCM TELEMETERED DATA 

by 

Dr. Paul Goodman 
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New Jersey 


PROBLEM STATEMENT 

The project objective was to create a method and equipment which would take as in- 
put a PCM signal containing data, and produce as output a single number (called K in the 
following) which would reliably predict the relative ease or difficulty of recovering error- 
free data from the tape. 

Such a method would allow greater care to be used when processing a tape predicted 
(by K) to be difficult. Conversely, poor data recovery from a tape predicted to be excel- 
lent would indicate malfunctioning processing (data recovery) equipment. 

Procedures and Methods 

A perfect, noiseless, split phase PCM signal has the following useful property: the 
time interval between successive negative to positive (or positive to negative) transitions 
is restricted to three values fixed by the bit transmission rate. For example, when one 
thousand bits per second are transmitted, theory predicts measured transition times of 
1, 1.5, and 2 milliseconds. Let us define three narrow intervals centered around these 
three times, and let K be the fraction of all measurements falling within all three inter- 
vals. Then it seems clear that K = 1.0 for a perfect signal, and K should decrease as 
signal distortion increases, or equivalently as signal to noise ratio decreases. It also 
seems clear that increasing signal distortion causes more errors in reading data. 
Therefore a reasonable hypothesis is that a correlation exists between K and error rate, 
which will allow K to be used as indicated in the introduction. 

The work performed at GSFC this summer consisted of defining the problem, planning 
an experiment to test the hypothesis just mentioned, setting up equipment and carrying out 
the experiment, and interpreting the results, hi more detail, the following was done. 

1. Eleven tapes containing data as telemetered from the OGO S-49 satellite were 
chosen for the experiment. The tapes spanned a range from very poor to excel- 
lent data recovery. 

2. Equipment was set up and transition times were measured automatically and re- 
corded. Each tape was sampled randomly, and approximately one thousand sam- 
ples per tape were recorded. 

3. The measurements were transferred to punched cards, and a computer was then 
used to find K for each tape, using intervals of ±5% around each of the three 
noise-free transition times. 

4. The values of K were compared with a prior knowledge of data errors in order 
to establish the degree of correlation. 
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RESULTS AND CONCLUSIONS 


Out of eleven tapes, one could not be used because of experimental errors. Of the 
remainer, six showed good correlation, two were border line, and two were poor. 

In some cases a tape was run twice, and K was computed separately for each run. 
Differences in K indicated that larger samples were called for, not a surprising conclu- 
sion when one takes into account that one tape may contain forty million bits. Sample 
size restricted to one thousand per tape because each measurement had to be transferred 
by hand to a punched card, the resulting bottleneck making more samples impractical. 

The correlation seems as good as possible, considering the small sample size. 
Another unknown factor is errors caused by a malfunctioning process line. One tape in 
particular was run twice, K was the same both times and was high enough to indicate 
very few errors, yet relatively many errors had been noted when the tape was processed. 

The final conclusion is that the results are promising and justify further experiments. 
Equipment must be devised to automate the sampling procedure, especially the process of 
getting the samples into final form suitable for computer input. This will allow larger 
samples to be used and many more tapes to be analyzed. 

A report now being written by the author contains more complete discussions, data, 
conclusions, and recommendations for future investigations. This report forms the in- 
troduction to the more complete report. 


Data Selection 

The following tapes were selected by M. Pasternack as a representative cross- 
section of signal characteristics, based on a prior knowledge of slippage rate (SR) and 
bit error rate (BER). 


Station # 
Tape # 
SR 
BER 


Each tape was rated for SR and BER as H (horrible), P (poor), F (fair), E (excellent). H 
possible, future investigations should replace these with numerical values for SR and BER, 
in order to permit an actual computation of correlation. Tape #425 was never delivered, 
and tape #14, Station 20 was not used for technical reasons. 


Table 1 
Tape Ratings 


16 

16 

19 

19 

19 

20 

20 

20 

20 

20 

20 

20 

14 

213 

425 

450 

491 

14 

633 

878 

915 

1418 

1621 

1667 

P 

- 

H 

P 

E 

- 

E 

E 

H 

H 

F 

H 

H 

- 

P 

H 

P 

E 

E 

E 

P 

H 

F 

F 


Equipment 

Tapes were played on a Consolidated Electro- dynamics Tape Recorder, Model VR- 
2600. The signal was taken in all cases from channel 2 containing the so called "raw 
data." It was then, in most cases, passed through a low-pass filter in order to eliminate 
tape deck noise, whenever noise and signal were sufficiently separated in frequency. 
Tape deck noise was always above forty thousand cps, and signal frequency could be 
lowered by changing tape speed. The filter output drove a Hewlett Packard 52436 Elec- 
tronic Counter whose reading was printed by a Hewlett Packard 560A Digital Recorder. 
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The counter had a 10 megacycle clock, permitting resolution to 0.1 microsecond, which 
was more than adequate. 

The numbers thus obtained were punched on cards and then sorted to give the num- 
ber of times each measurement occured. This was to be plotted as a distribution func- 
tion by the 4020 plotter, but the programs were never successful. However, the distri- 
butions were printed as output, and this permitted K to be computed by hand. Of course 
this could be done much more easily by the computer, but it was not clear exactly what 
was necessary when the program was written. 

Results and Conclusions 

Table 2 gives a summary of measurements, from which our conclusions are drawn. 
Let us make a very simple rule and see how well SI is predicted. We assume that H 
corresponds to 0 < K < 0.25, P to 0.25 < K < 0.5, F to .50 < K < .75, and E to 0.75 < K 
< 1*0. In cases where two runs correspond to one tape, the average value of K is taken. 

Table 2 

Computed values of K 


’ 

Tape 

SR 

BER 

K 

f 

f 

Run # 

19-450 

P 

H 

.21 

40 

8 

1 

19-14 

p~“ 

H 

.39 

60 

16 

2 

19-425 

H 

P 

.27 

40 

8 

3 

20-633 

E 

E 

00 

40 

1 

4 

20-1667 

H 

F 

.61 

40 

1 

5 

20-878 

E 

E 

.86 

40 

8 

6 

20-915 

H 

P 

.85 

40 

8 

7 

20-915 

H 

P 

.85 

00 

8 

8 

20-1621 

F 

F 

.30 

00 

1 

9 

20-1621 

F 

F 

.50 

30 

1 

10 

19-491 

E 

P 

.74 

40 

8 

11 

19-991 

E 

P 

.34 

40 

8 

12 

20-1418 

H 

H 

.69 

40 

ar 

13 


* 
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Table 3 shows the results for this and for predicting SR by the same rule, 


Table 3 

Predicted and Actual BER, SR 


Tape 

K 

BER 

(Predicted) 

BER 

(Actual) 

SR 

(Predicted) 

SR 

(Actual) 

19-450 

.21 

H 

H 

H 

P 

16-14 

.39 

P 

H 

P 

P 

19-425 

.27 

P 

P 

P 

H 

20-633 

.81 

E 

E 

E 

E 

20-1667 

.61 

F 

F 

F 

H 

20-878 

.86 

E 

E 

E 

E 

20-915 

.85 

E 

P 

' E 

H 

20-1621 

.40 

P 

F 

P 

F 

19-491 j 

.54 

F 

P 

F 

E 

20-1418 

.69 

F 

H 

F 

R 


BER is predicted correctly for 5 tapes, in three cases the prediction must be 
shifted up or down by one category (i.e., from P to H for tape 16-14) and in two cases a 
shift of two categories is necessary. Corresponding numbers for SR are three correct, 
four must be shifted by one category two need a shift of two, and one prediction is as far 
as possible from being correct. Clearly there is a greater correlation between K and 
BER than K and SR. In the former case results seem quite promising, leading to a ten- 
tative conclusion that K can be used to predict BER. 

Returning to Table 2, we see that three tapes were run twice (runs 7-12). If K is to 
have any meaning it must be constant for any one tape. Therefore, the changes in K be- 
tween runs 9 and 10, and between 11 and 12 must be eliminated. The cause of such 
changes is very likely the very limited sample size, and the fact that each set of samples 
were taken from, at most, about 25% of the tape. The number of bits per tapes may be of 
the order of ten million, thus, one thousand samples is not very much in comparison. 
Therefore, future experiments must use many more samples, and the samples should be 
spread uniformly over the entire part of the tape which contains data. 

Visual Inspection of Waveforms 

A visual inspection of signals on all eleven tapes was made. The object was to ex- 
plore the possibility of direct prediction of SR or BER by some easily recognized pattern 
or signal distortion. The tapes were divided into four categories on die basis of observed 
signal to noise ratio, and the results are shown in Table 4. 

Comparing this with Table 3 we see that BER was predicted correctly five times by 
K vs. four times visually, SR was correct three times in either case. Apparently K still 
is more indicative of BER than SR, and is a better predictor than visual recognition. 

Other reasons for preferring K are: it does not depend on human interpretation, it may 
lead to more sophisticated ways to treat the samples in a statistical sense. 
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Table 4 

Visual Rating of Signals 


Tape 

Rating 

(Visual) 

BER 

SR 

19-450 

H 

H 

H 

16-14 

P 

P 

H 

19-25 

F 

P 

P 

20-633 

E 

E 

E 

20-1667 

E 

F 

F 

20-878 

P 

E 

E 

20-915 

P 

E 

P 

20-1621 

E 

P 

F 

19-491 

F 

F 

P 

20-1418 

P 

F 

H 


The waveforms were reduced to strip charts showing amplitude vs. time, using two 
tape recorders, a Consolidated Electrodynamics VR-2600, and an Ampex FR-600. Each 
tape was played on the VR 2600 at 3.75 i.p.s. and recorded on the FR 600 at 60 i.p.s., then 
the recording was played back at 1-7/8 i.p.s. and recorded on a Honeywell 1508 Visi- 
corder. ft was necessary to use FM modules on the tape recorders for good low-frequency 
response. 

Future Work 

Before this work is continued on an experimental basis, it is necessary to set up 
equipment capable of automatically transferring readings from the HP electronic counter 
to either punched cards or tape, the latter being preferable because of the rapidity with 
which it can be entered into a computer. A larger number of tapes should be chosen, and 
numerical values of BER should be known for them. Then if K and numerical BER are 
related, it will be apparent from experimental results. A range of ±5% for computing K 
is not necessarily best, and a range of such values should be used to see the effect of 
change. 

Theoretical investigations should be made for comparison with experimental results. 
Two questions of a statistical nature seem appropriate: 

1. What is the expected distribution of transition times as a function of signal to 
noise ratio, for various well known types of noise such as white noise, shot 
noise, etc? 

2. What is an appropriate number of samples to ensure K being statistically 
reliable? 

These are questions the author hopes to investigate during the coming school year. 
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DEVELOPMENT OF THE 10.6-MICRON LASER 

by 

Dr. Zbigniew D. Jastrzebski, Professor 
Lafayette College, Easton, Pa. 


STATEMENT OF THE PROBLEM 

The use of 10.6-micron laser for optical communication systems may offer certain 
advantages as compared with the pulsed ruby laser used presently for tracking experi- 
ments. These advantages are higher degree of beam coherence, less energy required 
for an information bit at the receiver, and less atmospheric loss per unit distance. 

The GSFC Optical Systems Branch has initiated the program aiming at the develop- 
ment of an effective optical tracking system for communication between the Goddard 
Optical Research Facility and the Echo n satellite. As a part of this program, I was 
assigned to study and develop the 10.6-micron nitrogen-carbon dioxide laser. 


PRESENT STATE OF KNOWLEDGE 

A strong C. W. lasing action in a nitrogen-carbon dioxide mixture under low pres- 
sures (''1mm Hg) has been reported in literature. Lasing in infrared takes place when 
vibrationally excite N 2 molecules (V = 1) excite through inelastic collisions the C0 2 
molecules (00°0) to the upper laser level (00°1). Nitrogen molecules fall back to their 
vibrational ground-level (V = 1), leaving behind a large population of the excited C0 2 
molecules. The main output lines occur at 10.5915 microns and at 10.5716 microns on 
two rotational transitions of the 00° 1 - 10°0 vibrational band with a C.W. power output 
which may be as high as 11.9 watts. 


EXPERIMENTAL 

Two experimental setups have been designed which differ mainly in the arrangement 
of flat reflective mirrors in respect to the Vycor discharge tube of 24mm I.D. and 122 cm 
length. In one setup the goldplated mirrors form an integral part of the laser assembly, 
being vacuum-tight connected to the ends of the discharge tube through metallic bellows. 
In the second setup the gold-plated mirrors are located externally to the discharge tube 
which is terminated by vacuum-tight Brewster angle windows. These windows and one 
of the two gold-coated mirrors are made of Barium Fluoride. This mirror has one- 
millimeter diameter hole in the gold coating for coupling out the radiation; the other 
mirror is 100% opaque. 

The discharge tube is connected through an outlet with the manifold of the vacuum 
system so that it is possible to control precisely low partial pressures of the gaseous 
ingredients. The necessary pumping power is obtained using a Viking Invader Trans- 
mitter (C.W. 1000 watts) providing rf discharge of 28 Mc/sec. Only qualitative meas- 
ures for the detection of the infrared laser beam are now available. 
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CONCLUSIONS 

The experimental setups have been found to function satisfactorily. At the moment 
of this writing, no appreciable lasing action has yet been noticed, but concrete results 
are expected very shortly. 

This project has been carried out under the guidance of Mr. Nelson McAvoy and 
with assistance of the other staff members of the Quantum Optics Section. 
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